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Abstract — Several models for the simulation of photon elastic 
scattering are quantitatively evaluated with respect to a large 
collection of experimental data retrieved from the literature. 
They include models based on the form factor approximation, 
on S-matrix calculations and on analytical parameterizations; 
they exploit publicly available data libraries and tabulations 
of theoretical calculations. Some of these models are currently 
implemented in general purpose Monte Carlo systems; some have 
been implemented and evaluated for the first time in this paper 
for possible use in Monte Carlo particle transport. The analysis 
mainly concerns the energy range between 5 keV and a few MeV. 
The validation process identifies the newly implemented model 
based on second order S-matrix calculations as the one best 
reproducing experimental measurements. The validation results 
show that, along with Rayleigh scattering, additional processes, 
not yet implemented in Geant4 nor in other major Monte Carlo 
systems, should be taken into account to realistically describe 
photon elastic scattering with matter above 1 MeV. Evaluations 
of the computational performance of the various simulation 
algorithms are reported along with the analysis of their physics 
capabilities. 

Index Terms — Monte Carlo, simulation, Geant4, X-rays 



I. Introduction 

PHOTON elastic scattering is important in various exper- 
imental domains, such as material analysis applications, 
medical diagnostics and imaging [1]; more in general, elastic 
interactions contribute to the determination of photon mass 
attenuation coefficients, which are widely used parameters in 
medical physics and radiation protection [2]. In the energy 
range between approximately 1 keV and few MeV, which is 
the object of this paper, the resolution of modern detectors, 
high intensity synchrotron radiation sources and, in recent 
years, the availability of resources for large scale numeri- 
cal computations have concurred to build a wide body of 
knowledge on photon elastic scattering. Extensive reviews of 
photon elastic scattering, that cover both its theoretical and 
experimental aspects, can be found in the literature (e.g. [3]- 
[5]). 

This paper addresses this topic under a pragmatic perspec- 
tive: its simulation in general purpose Monte Carlo codes 
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for particle transport. Photon interactions with matter, both 
elastic and inelastic, play a critical role in these systems; their 
modeling presents some peculiarities, because the software 
must satisfy concurrent requirements of physical accuracy and 
computational performance. 

Photon-atom elastic scattering encompasses various interac- 
tions, but Rayleigh scattering, i.e. scattering from bound elec- 
trons, is the dominant contribution in the low energy regime 
and, as energy increases, remains dominant in a progressively 
smaller angular range of forward scattering. Rayleigh scat- 
tering models are implemented in all general-purpose Monte 
Carlo systems; comparison studies have highlighted discrepan- 
cies among some of them [6], nevertheless a comprehensive, 
quantitative appraisal of their validity is not yet documented 
in the literature. It is worthwhile to note that the validation of 
simulation models implies their comparison with experimental 
measurements [7] ; comparisons with tabulations of theoretical 
calculations or analytical parameterizations, such as those that 
are reported in [8] as validation of Geant4 [9], [10] photon 
cross sections, do not constitute a validation of the simulation 
software. 

This paper evaluates the models adopted by general-purpose 
Monte Carlo systems and other modeling approaches not yet 
implemented in these codes, to identify the state-of-the-art 
in the simulation of photon elastic scattering. Computational 
performance imposes constraints on the complexity of physics 
calculations to be performed in the course of simulation: 
hence the analysis is limited to theoretical models for which 
tabulations of pre-calculated values are available, or that are 
expressed by means of simple analytical formulations. To be 
relevant for general purpose Monte Carlo systems, tabulated 
data should cover the whole periodic table of elements and 
an extended energy range. The accuracy of elastic scatter- 
ing simulation models is quantified with statistical methods 
comparing them with a wide collection of experimental data 
retrieved from the literature; the evaluation of physics capabil- 
ities is complemented by the estimate of their computational 
performance. These results provide guidance for the selection 
of physics models in simulation applications in response to 
the requirements of physics accuracy and computational speed 
pertinent to different experimental scenarios. 

Special emphasis is devoted to the validation and possible 
improvement of photon elastic scattering modeling in Geant4; 
nevertheless, the results documented in this paper can be of 
more general interest also for other Monte Carlo systems. 
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II. Physics overview 

Photon-atom elastic scattering has been the object of the- 
oretical and experimental interest for several decades; only a 
brief overview is summarized here to facilitate the comprehen- 
sion of the software developments and simulation validation 
results documented in this paper. 

Conventionally, photon-atom interactions are classified into 
elastic and inelastic processes, depending on whether the 
photon energy changes in the center of mass frame: although 
commonly applied in atomic physics practice and Monte Carlo 
transport, this distinction is somewhat arbitrary, since QED 
(quantum electrodynamics) radiative corrections and target 
recoil actually make all processes inelastic and, on the experi- 
mental side, the classification of a process as elastic or inelastic 
is limited by source bandwidth and detector resolution. 

The physics models considered in this study concern iso- 
lated photon-atom scattering: they neglect the effects of the 
environment of the target atom and multi-photon effects as- 
sociated with incident beams. They do not take into account 
either any internal degrees of freedom of the initial or final 
atom, such as its orientation; therefore they do not describe 
X-ray scattering from electrons in molecules and oriented 
solids. This scenario is consistent with the assumptions of 
general-purpose Monte Carlo systems for particle transport, 
where particles are assumed to interact with free atoms of the 
medium. 

A. Theoretical calculations 

Elastic photon-atom scattering encompasses various types 
of interactions, which are usually distinguished as: 

• Rayleigh scattering, for scattering from the atomic elec- 
trons, 

• nuclear Thomson scattering, for scattering from the nu- 
cleus considered as a point charge, 

• Delbriick scattering, for scattering from the field of the 
nucleus, 

• nuclear resonance scattering, for scattering from the in- 
ternal structure of the nucleus. 

These amplitudes are coherent, i.e. not physically distin- 
guishable; the total elastic amplitude is commonly calculated 
as the sum of these separate amplitudes. There are still 
unresolved questions about the additivity of the amplitudes 
from composite systems, due to the complex interplay in a 
composite system between partitioning the total amplitude and 
the reduction from a many -body photon-atom interaction to an 
effective single particle formalism [11]; open issues concern 
the effects of electron correlations and of the finite lifetimes 
of atomic excited states [3]. 

Rayleigh scattering calculations are commonly performed 
by means of an approximation, that describes the scattering 
amplitude in terms of the form factor, which represents the 
Fourier transform of the charge density of the atom. This can 
be simply understood considering the two Thomson ampli- 
tudes from a free (point-like) particle [12]: 

Af { = -r cos , A^ = -r (1) 



leading to a differential cross section 

^ = =f<lW«) (2) 

which, for coherent scattering, is modified in presence of a 
charge distribution in 

^ = r ^(l + cos^)F\ q ,Z) (3) 

where a represents the cross section, Q is the solid angle, ro is 
the classical electron radius, 6 is the scattering angle, F(q, Z) 
is the atomic form factor as a function of momentum transfer q 
and atomic number Z of the target atom. This approximation is 
valid for photon energies much greater than electron binding 
energies and for non-relativistic momentum transfers, which 
in this case is meant as hq <C m e c, where H is the reduced 
Planck constant, m e is the electron mass and c is the speed of 
light. In the so-called "modified form factor" approximation 
each subshell charge density is multiplied by a different factor. 

Various calculations of non-relativistic, relativistic and mod- 
ified relativistic form factors are documented in the literature; 
the tabulations by Hubbell et al. [13], Hubbell and 0verb0 
[14], and Schaupp et al. [15], respectively of non-relativistic, 
relativistic and modified form factors, are representative of this 
approach. 

For forward scattering the difference between the form 
factor amplitude and the exact amplitude is described by two 
anomalous scattering factors [3], that take into account the 
persistence of binding effects even in the high energy limit of 
Rayleigh scattering. 

A more consistent approach for the description of elastic 
photon-atom scattering involves the evaluation of the relativis- 
tic second order S -matrix element in independent particle ap- 
proximation. Relativistic quantum electrodynamics, treated in 
lowest non-vanishing order in e 2 , provides the basic theoretical 
framework for these calculations; neglect of higher order terms 
in e 2 in the calculation means that radiative corrections are not 
taken into account. The mass of the atomic nucleus is assumed 
to be infinite compared to the photon energy in question; 
therefore only scattering from atomic electrons is considered. 
This many-electron state is treated in a self-consistent central 
field approximation, in which each electron is assumed to 
move in an average potential due to all atomic electrons and 
the nucleus. 

Numerical evaluation of the second order S -matrix for single 
electron transitions in a potential was first attempted in the 
1950s by Brown and co-workers [16]. This calculation scheme 
requires considerable computing resources, since the multipole 
expansion of the photon field converges slower and slower 
for increasing energies; systematic evaluation of the S-matrix 
element has become practical relatively recently, thanks to 
wide availability of large scale computing facilities. Extensive 
calculations based on the S-matrix approach (e.g. [11], [17], 
[18]) have been performed by Kissel, Pratt and co-workers. 

B. Data libraries 

Some results of theoretical calculations of photon elastic 
scattering are publicly available in the form of tabulations 
distributed as data libraries. 
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EPDL97 [19] (EPDL Evaluated Photon Data Library, 1997 
version), which is part of the ENDF/B-VII.l [20] evaluated 
nuclear data file, includes total cross sections, form factors and 
anomalous scattering factors for atoms with atomic number 
between 1 and 100, and for photon energies from 1 eV to 
100 GeV. The form factors in EPDL97 are the non-relativistic 
ones calculated by Hubbell [13]; the anomalous scattering 
factors are by Cullen [21]; the total cross sections derive 
from calculations combining Thomson scattering, form factors 
and anomalous scattering factors, which were numerically 
integrated. 

EPDL97 documentation reports rough estimates and quali- 
tative comments about the accuracy of the tabulated data, but 
it does not document how these estimates were produced. To 
the best of our knowledge systematic, quantitative validation 
of EPDL97 coherent scattering data is not documented in the 
literature. 

EPDL97 is extensively used in Monte Carlo simulation; 
details are given in section |rQ| 

The RTAB [22] database encompasses a set of tabulations of 
photon elastic scattering cross sections, which are the result of 
various methods of calculation, and of components for their 
calculations, such as form factors and anomalous scattering 
factors. The differential cross sections listed in RTAB are 
based on: 

• numerical S -matrix calculations by Kissel and Pratt, 

• relativistic form factors, 

• non-relativistic form factors by Hubbell et al. [13], 

• modified relativistic form factors, 

• modified relativistic form factors with angle-independent 
anomalous scattering factors. 

Apart from the data derived from other sources, such as 
Hubbell 's form factors, all the data in the RTAB database have 
been consistently computed in the same Dirac-Slater potential. 
The tabulations have been generated on a 97-point grid for 
scattering angles between and 180 degrees, and on a 56- 
point grid for energies between 54.3 eV and 2.7541 MeV; 
they are listed for atomic numbers from 1 to 100. 

Two sets of cross sections based on S-matrix tabulations 
are available in RTAB: one takes into account both the 
Rayleigh and Thomson scattering amplitudes, which are added 
coherently, and one is limited to the Rayleigh scattering am- 
plitude. In the following, "S-matrix" or "SM" identify RTAB 
tabulations involving both amplitudes; whenever tabulations 
limited to Rayleigh scattering are used, they are explicitly 
identified in the text. 

Comparisons of measured data with S-matrix calculations 
are reported in some experimental publications, which are lim- 
ited to the configuration (photon energy, scattering angle and 
target material) of the experiment they describe. A systematic, 
quantitative validation of RTAB data is not documented yet 
in the literature. This database has not been exploited yet in 
general purpose Monte Carlo systems. 

III. Photon elastic scattering in Monte Carlo 
codes 

General purpose Monte Carlo codes for particle transport 
include algorithms for the simulation of Rayleigh scattering, 



most of which are based on the form factor approximation; 
they do not appear to account for other amplitudes involved 
in photon elastic scattering, apart from the implicit inclusion 
of Thomson scattering in the calculation of total cross sections 
by the codes that use EPDL97. 

These codes assume that photons interact with free atoms 
of the medium: this assumption neglects that the molecular 
structure and the structure of the medium can affect coherent 
scattering. Some codes provide users the option to directly 
input tabulated Rayleigh scattering cross sections and form 
factors, which may account for molecular effects specific to 
a given material; such tabulations are available for a limited 
number of materials (e.g. [23]). Photon elastic scattering 
by molecules is not treated in this paper. The information 
summarized here concerns models implemented in Monte 
Carlo codes to describe elastic scattering by non-polarized 
photons; some Monte Carlo systems also include algorithms 
for Rayleigh scattering by polarized photons, which are not 
considered here. 

EGS5 [24] and EGS4 [25] calculate total Rayleigh scatter- 
ing cross sections from the tabulations by Storm and Israel 
[26], which in turn derive from the integration over angle 
of equation ([2]), where form factors by Hanson [27] were 
used. They sample the coherent scattering angle based on the 
relativistic form factors by Hubbell and 0verb0 [14]. EGSnrc 
[28] can use three sets of cross section data: those by Storm 
and Israel as EGS4, and the additional options of EPDL97 and 
XCOM [29] cross sections. 

ETRAN [30] uses the cross sections of the XCOM database, 
which are based on the relativistic form factors of Hubbell 
and 0verb0 [14], and samples the change in photon direction 
according to the form factor approximation. 

FLUKA [31], [32] simulates Rayleigh scattering based on 
EPDL97. 

ITS [33] calculates total Rayleigh scattering cross sections 
based on the relativistic form factors by Hubbell and 0verb0 
[14], but it samples the scattering angle based on the non- 
relativistic form factors by Hubbell et al. [13]. 

MCNP5 [34] and MCNPX [35] neglect Rayleigh scattering 
when the "simple option" of photon transport is chosen; 
otherwise they also base the simulation on the form factor 
approximation. They provide different options of Rayleigh 
scattering data, the most recent of which uses EPDL97; the 
other options use the older EPDL89 [36] version of EPDL, 
data from ENDF/B-IV [37] based on non-relativistic form 
factors by Hubbell et al. [13], and the tabulations by Storm 
and Israel [26], the latter limited to a few elements with 
atomic number greater than 83. Nuclear resonance fluores- 
cence, nuclear Thomson scattering and Delbriick scattering are 
not treated [38]. 

Penelope [39] 2008 [40] and 2011 [41] versions calculate 
total Rayleigh scattering cross sections and scattering angles 
based on EPDL97 [19]; earlier versions exploited analytical 
approximations [42] to Hubbell 's non-relativistic form factors 
for the calculation of differential and total cross sections. 
While Penelope 2008 and 2011 documentation [40], [41] 
states that the total cross sections and form factors used in 
Penelope are from EPDL97, the tabulations of these quantities 
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distributed with the Penelope code appear different from those 
in EPDL97: presumably, the tabulations included in Penelope 
have been recalculated, or interpolated, over a different energy 
or momentum transfer grid. 

GEANT 3 [43] handled Rayleigh scattering according to 
empirical formulae derived from EGS3 [44]; they consist of 
polynomial fits to cross sections by Storm and Israel, and 
relativistic form factors by Hubbell and 0verb0 [14]. 

The Geant4 toolkit encompasses various implementations of 
Rayleigh scattering. The latest versions at the time of writing 
this paper are Geant4 9.4p03 and Geant4 9.5, both released in 
December 2011 (9.4p03 being the most recent). The Geant4 
G40pRayleigh class implements a Rayleigh scattering model 
only applicable to a particle type identified in Geant4 as 
"optical photon", which is an object of the G40pticalPhoton 
class; this process is not considered in the following analysis, 
that deals with Geant4 Rayleigh scattering models concerning 
ordinary photons, which are instances of the G4Gamma class. 

Functionality for the simulation of Rayleigh scattering of 
ordinary photons was first introduced [45] in the low energy 
electromagnetic package [46], [47] of Geant4 0.1 version, 
based on total cross sections and non-relativistic form factors 
tabulated in the EPDL97 [19] data library; this modeling 
approach is identified in Geant4 9.5 as "Livermore Rayleigh 
model". 

Geant4 also includes two implementations of Rayleigh scat- 
tering reengineered from the Penelope [39] Monte Carlo code, 
respectively equivalent to the Rayleigh scattering algorithms 
in Penelope 2001 [48] and 2008 [40] versions, the latter also 
based on EPDL97 as the original Rayleigh scattering im- 
plementation available in Geant4. The reengineered Penelope 
2008 code uses the tabulations of Rayleigh scattering total 
cross sections and form factors distributed with Penelope, 
which, as discussed above, appear different from the native 
EPDL97 values used by the "Livermore Rayleigh" implemen- 
tation. 

A further Rayleigh scattering simulation model, imple- 
mented in the G4XrayRayleighModel class, has been first 
released in Geant4 9.5 version: it is defined in the software 
release notes as "implementing simplified Rayleigh scatter- 
ing", but no further information could be found in Geant4 
documentation, nor in the literature, about the physics model 
underlying this implementation. From the C++ source code, 
the authors of this paper could evince that this model cal- 
culates the total Rayleigh scattering cross section according 
to an analytical formula involving numerical parameters, and 
samples the scattering angle according to the distribution of 
Thomson scattering from a point-like charge. 

IV. Strategy of this study 

An extensive set of simulation models, which are repre- 
sentative of the variety of theoretical approaches documented 
in the literature, have been evaluated to identify the state-of- 
the-art of photon elastic scattering in the context of Monte 
Carlo particle transport. The physics models considered in this 
analysis involve the implementation of simple formulae in the 
simulation software or exploit available tabulations of complex 
theoretical calculations. 



The models for photon elastic scattering simulation evalu- 
ated in this paper concern total and differential cross sections: 
in particle transport, the former are relevant to determine the 
occurrence of the scattering process, while the latter determine 
the actual outcome of the scattering by defining the angular 
distribution of the scattered photon at a given energy and for 
a given target. 

All the models subject to study have been implemented in a 
consistent software design, compatible with the Geant4 toolkit, 
which minimizes external dependencies to ensure the unbiased 
appraisal of their intrinsic capabilities. 

A wide set of experimental data of photon elastic scattering 
has been collected from the literature for this study; simulation 
models are validated through comparison with these measure- 
ments. The compatibility with experiment for each model, and 
the differences in compatibility with experiment across the 
various models, are quantified by means of statistical methods. 

The measurements of photon elastic scattering reported 
in the literature mainly concern differential cross sections; 
therefore the validation process is focused on this observable. 
Due to the scarcity of total cross section measurements in the 
literature, the capability of the simulation models to calculate 
total cross sections consistent with experiment can be directly 
tested only over a small data sample; nevertheless, it can be 
indirectly inferred for total cross sections deriving from the 
integration of validated differential cross sections. 

The validation of the physics capabilities of the simulation 
models is complemented by the evaluation of their computa- 
tional performance. 

A. Simulation models 

The simulation models of differential cross section analyzed 
in this study implement calculations based on second order 
S -matrix and on the form factor approximation; they are 
summarized in Table |l] Some of these simulation models 
represent novel approaches with respect to those so far avail- 
able in Geant4 and in other general purpose Monte Carlo 
codes; among them, the paper examines whether a model 
based on S -matrix calculations, which constitutes a more 
consistent theoretical approach, but requires more complex 
computational operations, would be sustainable in the context 
of Monte Carlo particle transport. 

Various simulation models are based on the non-relativistic 
form factors by Hubbell et al. [13]: they are identified in 
Table |T| as "EPDL", "NF", "Penelope 2001" and "Penelope 
2008". They use different tabulations of the form factors 
as a function of momentum transfer, or different algorithms 
to sample the scattered photon direction based on the form 
factors. The model identified as "EPDL" uses non relativistic 
form factors as they are tabulated in EPDL97; it corresponds to 
the original Rayleigh scattering simulation method in Geant4 
low energy electromagnetic package [45], which has been 
reengineered in the context of the policy-based class design 



described in section IV-B The models identified in this paper 
as "Penelope 2001" and "Penelope 2008" correspond to al- 
gorithms reengineered from the respective Penelope versions: 
as mentioned in section [In| the Penelope 2008 model is 
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TABLE I 

Simulation models of photon elastic scattering: differential cross sections 



Identifier 


Underlying theoretical approach 


Tabulated data source 


EPDL 

EPDLASF 

Penelope 2001 

Penelope 2008 

RF 

NF 

MF 

RFASF 
MFASF 
SM 


Non-relativistic form factors 

Non-relativistic form factors + anomalous scattering factors 
Parameterization of non-relativistic form factors (as in Penelope 2001) 
Non-relativistic form factors (as in Penelope 2008-2011) 
Relativistic form factors 
Non-relativistic form factors 
Modified form factors 

Relativistic form factors and anomalous scattering factors 
Modified form factors and anomalous scattering factors 
Second order S-matrix 


EPDL97 
EPDL97 

Retabulated EPDL97 
RTAB 

Hubbell 1975, as reported in RTAB 
RTAB 
RTAB 
RTAB 
RTAB 



based on the EPDL97 data library, although re-tabulated over 
a different grid, while the Penelope 2001 model exploits 
analytical approximations to Hubbell 's non-relativistic form 
factors. The model identified in this paper as "NF" exploits 
the non-relativistic form factors listed in RTAB. 

Models based on relativistic and modified form factors, 
based on RTAB tabulations, are respectively identified as "RF" 
and "MF". To the best of the authors' knowledge, modified 
form factors have not been used yet in general purpose Monte 
Carlo codes. 

Additionally, calculations accounting for anomalous scatter- 
ing factors along with form factors are evaluated: the "RFASF" 
and "MFASF" models are based on RTAB tabulations, while 
the "EPDLASF" model exploits the anomalous scattering 
factors included in EPDL97, which have not been used yet 
in any of the general-purpose Monte Carlo systems relying on 
EPDL97. 

The model identified as "SM" is based on RTAB tabulations 
of second order S-matrix calculations. This model accounts for 
the sum of two coherent amplitudes: the so-called Rayleigh 
amplitude and nuclear Thomson scattering. A model based 
on S-matrix calculations limited to the Rayleigh scattering 
amplitude is discussed in section |V-D| To the best of the 
authors' knowledge, this paper documents the first implemen- 
tation based on S-matrix calculations for the simulation of 
photon elastic scattering by a general-purpose Monte Carlo 
system. 

Total cross sections are calculated corresponding to all the 
models listed in Table [TJ those that are highlighted as most 
relevant by the validation analysis of differential cross sections 
in section IV-BI are discussed in detail in section IV-CI Further- 
more, the total cross sections calculated by interpolation of 
Storm and Israel [26] tables, by the XCOM Photon Cross 
Section Database distributed by NIST (National Institute of 
Standards), by Brennan and Cowan [49] parameterizations of 
McMaster's [50] cross sections (commonly used in photon 
science), and by Geant4 G4XrayRayleighModel are evaluated. 
The energy ranges covered by these compilations are respec- 
tively 1-100 MeV (Storm and Israel), 0.03-695 keV (Brennan 
and Cowan) and 1 keV to 100 GeV (XCOM). The total cross 
section models examined in this paper are summarized in 
Table El 

Tabulations of total Rayleigh scattering cross sections are 
directly available in EPDL97, XCOM and [26]; tabulations for 
other physics models have been produced on the same energy 



grid as EPDL97 to facilitate the comparison of the various 
models. The tabulations associated with the simulation models 
based on RTAB listed in Table [I] were produced by integrating 
the corresponding differential cross sections over the angle. 



B. Software environment 

All the physics models evaluated in this paper have been 
implemented in the same software environment, which is 
compatible with Geant4; computational features specific to 
the original physics algorithms have been preserved as much 
as possible. The uniform software configuration ensures an 
unbiased appraisal of the intrinsic characteristics of the var- 
ious physics models. The correctness of implementation has 
been verified prior to the validation process to ensure that 
the software reproduces the physical features of each model 
consistently. 

The software adopts a policy -based class design [51]; this 
technique was first introduced in a general-purpose Monte 
Carlo system in [52]. This programming paradigm supports the 
provision of a series of alternative physics models for Monte 
Carlo transport, characterized by high granularity, which can 
be used interchangeably with great versatility, without im- 
posing the burden of inheritance from a pre-defined inter- 
face, since policies are syntax-oriented, rather than signature- 
oriented. 

Two policies have been defined for the simulation of 
photon elastic scattering, corresponding to the calculation of 
total cross section and to the generation of the scattered 
photon; they conform to the prototype design described in 
[54], [55]. A photon elastic scattering process, derived from 
the G4VDiscreteProcess class of Geant4 kernel, acts as a 
host class for these policy classes. All the simulation mod- 
els implemented according to this policy-based class design 
are compatible for use with Geant4, since Geant4 tracking 
handles all discrete processes polymorphically through the 
G4VDiscreteProcess base class interface. 

A single policy class calculates total cross sections for all 
the physics models that exploit tabulations; alternative tabula- 
tions, corresponding to different physics models, are managed 
through the file system. Specific policy classes implement the 
analytical calculations of total cross sections as in Penelope 
2001 and G4XrayRayleighModel. 

Two alternative policy classes sample the direction of the 
scattered photon utilizing the form factor approximation: they 
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TABLE II 

Simulation models of photon elastic scattering: total cross sections 



Identifier 


Method 


Source 


EPDL 

Penelope 2001 

MFASF 

SM 

XCOM 
Storm 
Brennan 
G4std 


Non-relativistic form factors with anomalous scattering factors 
Analytical model 

Integration of differential cross sections based on modified form factors with anomalous scattering factors 
Integration of cross sections based on second order S-matrix 
Based on form factor approximation 
Storm and Israel tabulations 

Brennan and Cowan's parameterization of McMaster et al. cross sections 
Parameterization in the Geant4 G4XrayRayleighModel class 


EPDL97 
Penelope 2001 
RTAB 
RTAB 
NIST 
[26] 
[49] 

Not documented 



differ by the sampling algorithm, respectively based on the in- 
verse transform method and on acceptance-rejection sampling 
[53]. Both can be used with any tabulated form factors, which 
are managed through the file system. 

A dedicated policy class deals with the generation of the 
scattered photon based on S-matrix calculations; the scattering 
angle is sampled based on the inverse transform sampling 
method. 

Dedicated policy classes generate the scattered photon ac- 
cording to the analytical formulations used in Penelope 2001 
and G4XrayRayleighModel. 

The final state generation based on the form factor approxi- 
mation involves one-dimensional interpolations only, while the 
use of S-matrix calculations requires bi-dimensional interpo- 
lation of RTAB data. Interpolation algorithms are discussed in 
section IV-FI 

The software design adopted in this study ensures greater 
flexibility than the design currently adopted in Geant4 electro- 
magnetic package, since it allows independent modeling (and 
test) of total cross section calculation and photon scattering 
generation. Since policy classes are characterized by a single 
responsibility and have minimal dependencies on other parts of 
the software, the programming paradigm adopted in this study 
also facilitates the validation of physics models: validation 
tests involve only the instantiation of an object of the policy 
class and the invocation of the member function implementing 
the policy, therefore they expose only the intrinsic features 
of the physics models, excluding possible effects on physics 
behaviour due to other parts of the Monte Carlo system. 
Equivalent validation tests of physics models as currently 
implemented in Geant4 require instead a full scale simulation 
application, including the creation of a geometry model even 
when it is conceptually redundant (e.g. to validate a cross 
section calculation). 

C. Experimental data 

Experimental data for the validation of the simulation mod- 
els were collected from a survey of the literature. 

The sample of experimental differential cross sections con- 
sists of approximately 4500 measurements [3], [58]- [158], 
which concern 69 target atoms and span energies from 
5.41 keV to 39 MeV, and scattering angles from 0.5 to 165 
degrees. An overview of this data sample is summarized in 
Tables Hill and El 

Experimental total cross section data are scarce in the litera- 
ture. To the best of our efforts we could retrieve only a journal 



publication [159] reporting a measurement at 661.6 keV for 
four elements: barium, tungsten, lead and bismuth. A larger 
set of measurements by the same group as [159], concerning 
18 target elements and two photon energies (279.2 keV and 
661.6 keV), is documented in [160]; this reference includes 
the four values published in [159]. The experimental errors 
reported in [160] vary between 3% and 6% for the mea- 
surements at 661.6 keV, and between 0.3% and 0.5% for 
the data at 279.2 keV: since the reported precision of the 
two sets of measurements in similar experimental conditions 
differs by an order of magnitude, and better than 1 % precision 
appears inconsistent with typical experimental errors in the 
field documented in the literature, one may wonder whether 
some of the experimental uncertainties listed in [160] could 
be affected by a typographical slip. 

It is worthwhile to note that the available experimental 
differential cross section data are too scarce to calculate total 
cross sections based on them. Therefore they could not be 
exploited for the validation of theoretical total cross section 
models. 

A small sample of experimental measurements of photon 
elastic scattering at the electronvolt scale has been retrieved; 
data at such low energy are scarce in the literature, and consist 
mostly of measurements of molecular cross sections, which 
are outside the scope of this paper. A few measurements 
of absolute total cross sections, that concern atomic gaseous 
targets (mostly noble gases) in the energy range of approx- 
imately 1-10 eV, have been retrieved in the literature: the 
intrinsic physical characteristics of these data are congruous 
for comparison with calculations performed in independent 
particle approximation, since molecular effects, which become 
important at low energies, are expected to be minimized in 
atomic gaseous targets. Some of these data involved normal- 
izations to other values to obtain absolute cross sections: these 
manipulations can be source of systematic effects and, if the 
normalization is based on theoretical references, would make 
the data inappropriate for the validation of simulation models. 
Although this very low energy data sample is too limited to 
serve for proper simulation validation, it enables at least a 
qualitative appraisal of EPDL97 cross sections, that extend 
down to 1 eV. 

Some experimental cross sections have been published only 
in graphical form; numerical values were digitized from the 
plots by means of the Engauge [56] software. The error 
introduced by the digitization process was estimated by dig- 
itizing a few plots, whose numerical content is explicitly 
reported in the related publications. The reliability of the 
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TABLE III 



Summary of the experimental data used in the validation analysis: atomic numbers 3-72 



Element 


Ih y\ f* vcy\i rQnfrp 
X-/llClg,y IdllgjC 


A IT (1 1 f * fQIKTP 

/A.llg,lC lculg,C 


OcullJJlC 


rvcici Climes 


z 


vJ J LLLUKJl 


v Ke v; 


^U-CgjlCCS ^ 






3 


Li 


1600 


124 


I 


rs8i 

PoJ 


4 


DC 


SQ S4 


4S 
4J 


Q 

y 


n o^i 


S 




145 4 


9 0^ 


4 


TQ1I T147I 






99 1 1 /l^ ZL 
ZZ. 1— 14J .4 


9 0^ 1 ^ 

Z.U3— 1 33 


1 1 
i j 


ron rQ7i n 1 n n4n ri47i 


7 

/ 


N 


145 4 


9 0^ 

Z.U3 




n 47i 
Li4 / \ 


12 


Ma 


99 1 6 SQ S4 

ZZ. 10-J7.Jt 


QO 191 


4 


n^^i ri48i ri4Qi 

[1J>JJ, [14oJ, [14!7j 


13 


Al 




1.6-145 


86 


TS81 T661 T6Q1 T701 T761 T8^1 TQ11 TQ61 TQ81 nn^l T1071 T1971 n^91 n^Sl 
[joj, [uuj, [U"j, [/L/j, [ / uj , [ojj, [!?ij, [!?uj, ["oj, [il/jj, [iu/j, [iz/j, [ijzj — [ijjj, 

[139]-[142], [148], [149], [153], [156] 


14 


Si 


17 AA SQ S4 


90-121 


3 


TQ71 n^1l T14Q1 


16 




22.1-59.54 


121-133 


2 


[111], [149] 


Zl 


at 


1 7 44-^0 8S 


90-90 


2 


[1 D 1J, [1 JJJ 


99 
zz 


Ti 
11 


1 7 44 SQ S4 
1 / .44—9:7.94 


10 1 60 
1U— lOU 


1 8 
1 o 


r7^i ni^i ni^i ri97i n^n n^^i ri48i ri4Qi 

L/Jj, lLL3\, U*9J, llZ/\, LIjIJ, liJJ], li4o\, LL4^J 


9^ 

ZJ) 


v 

V 


17 44 SQ S4 


10 1 60 
iu— lOU 


1 7 

1 / 


r76i rQ7i ni^i ni^i n^n n^^i ri4Qi 


94 

Z4 




17 44 ^ 9Q 
1 / .44— DD.Ly 


QO QO 

yu— 


Q 


n^n n^^i n^si 

[191J, [i99J— [19JJ 


9S 

ZJ 


Mn 
1V1I1 


1 7 44 ^ 9Q 
1 / .44— DJ.Zy 


QO QO 


A 
u t 


n^n n^^i n^^i 

\_l^>l\, U991, 


96 

ZD 


re 


A 4 1 1 90 S 
0.4— 1 1ZU. J 


9 79 1^0 
Z. / Z-Iju 


^8 

JO 


T7S1 T811 T881 TQ71 T1 1 SI T1971 n^11 T1 ^41 M^SI ri411 T1481 T14Q1 
L/JJ, LolJ? L""J? \y ' h L^^9J, [1Z/J, [IjIJ, [1j4J, [IjjJ, [141J, [14oJ, [14yj 


97 
z / 


P n 
v_A) 


1 7 44 SQ S4 


10 1 60 
1U— lOU 


1 1 
19 


r 1 1 on ri97i ri^n n^si ri48i 

[119J, [1Z/J, [IjIJ, [IjjJ, [14oJ 


28 


Ni 


7 47 81 
/ .4 /— o 1 


4S 1 4S 
43—143 


99 
ZZ 


T661 T6Q1 T7S1 M 1 SI T1971 n^11 TL^I T1^S1 T1481 T14Q1 
[ooj, LO-^J? L'^J? LA 191, L-^^'J? LJ-91J, [IjjI, [IjJI, [i^oi, [i^^i 


9Q 




8 04 1 408 1 

O.U4 — 14UO. 1 


1 09 16S 

1 .UZ— 1DJ 


^06 


T64I T6SI T7^I T7SI T78I T81I T8^I T87I TQ1I TQ4I TQ6I TQ8I T10SI T10QI T114I 
[04j, LOJj, L'^J? L'^J? L'Oj, L"iJ, Lojj— [o/j, L-'iJ? L-^^j, l^oj— [^oj, [iujj, [iu^j, [ii4j, 

T11S1 T1991 T1941 T1971 T19Q1 n^11 n^SI n^QI T1491 T1471 T1 SOI T1 S^l T1 S61 
|_119J, [1ZZ1~ [1Z41, [IZ/l, [IZ^l, [1 j 11~ [1 J>J1, [1 Jy\— [14Z1, [14/1— [1JU1, L199J, [IJOl 


30 


Zn 


8 K\ 661 6 
O.kjD — UU1 .u 


10-165 


87 


T7S1 T771 T8^1 T8S1 T881 TQ^l TQ71 T1191 Tll^l T11S1 T1971 n^1l n^Sl T14S1 
L/jj, [/ /J, [ojj, L^JJ — L^oj, \y /j, [iizj, [iijj, [iijj, [iz/j, [ijij, [ijjj, [ih-jj, 

[148], [149], [152] 




Ha 


16.58 


90 




L1991 


^9 


UC 


Q 8 SQ S4 


QO 1 11 
y\J — Ijj 


4 


T7Si n 1 11 n 4Qi 


33 


As 


99 1 
ZZ. 1 


1 ^ 
199 


1 

1 


n 1 11 


34 


Se 


22.1 


133 


I 


n 1 1 1 
[ii 11 


36 


Kr 


59.54 


20 


11 


[151] 


37 


Rb 


22.1 


133 


I 


[111] 


38 


Sr 


14.14-30.85 


90-90 


10 


[75], [132], [134], [135] 


jy 


Y 


17 44 SQ S4 

1 / .44— J7. Jt 


10 1 60 


1 Q 

1 y 


ni^i n^n n^^i n^si n^Qi ri40i ri48i ri4Qi 

[11J>J, [1J>1J, [IJJJ, [IJJJ, [IJ^J, [14UJ, [14oJ, [14!7j 


40 


Zr 


1 ^ QS 661 6 

1 J.7J— DDI .O 


10 1 6S 


67 
o / 


T7SI T80I T8SI T87I TII^I T1 1 SI T1901 T19S1 T1971 T1 ^41 n^SI T1481 T14Q1 

[/JJ, L^UJ, [OJJ— [O/J, [llJj, [11JJ, [1ZUJ, [1ZJJ— [1Z/J, [1J>4J, [IJJJ, [14oJ, [14!7j 


41 


Nb 


13.95-59.54 


10-165 


66 


[75], [80], [85]-[88], [113], [115], [127], [131], [135], [148], [149] 


42 


Mo 


13.95-1408.1 


2-165 


220 


[64], [65], [67], [75], [76], [80], [85]-[87], [97], [98], [112]-[115], [121], [123], [124], 
n 97i ru^i n^si n40i n4Si n48i nsoi ns9i nssi ns6i 

[1Z/J, [iJJj, [IJJJ, [1H-L/J, [1H-JJ, [1H-OJ — [1JUJ, L A J-^J » L A J J J » L A JOJ 


44 

44 


IN. LI 


99 1 6 SQ S4 

ZZ. 1U — J". Jt 


90-121 


2 


T7S1 T14Q1 


4S 
4J 


Rh 


99 1 6 SQ S4 
ZZ. ID— 9:7.94 


117 191 
1 1 /— 1Z1 


9 
z 


T148I T14QI 


46 
40 




< 41 SQ S4 
J .4 l-J^Jt 


QO 1 ^0 
y\J— 1 3U 


1 8 
1 o 


n 1 si ri97i n^n n^^i n^si n^7i ri48i ri4Qi 

[119J, [IZ/l, [1J>11, [1 DD\— [1 JJ], [1J>/1, [14ol, [14^1 


Al 
4 / 


A (T 


9.41— OOl .o 


S 1 9 1 6S 
3 . 1 Z— 1 09 


1 9^ 

1ZJ 


TSQI T60I T6^I T6SI T7SI T8SI T87I TQ4I T1 1 11 ni^I ri9^I T194I n^1l n^^I 
Ljyj, LOUJ, L09J-L09J, l/D], Lojj-Lo/J, \y4\, L a 1 1J, LH9J, lLZ3\, LLZ4J, IL3L\, iL33\- 

n^si n^7i ri4Qi nsoi 

[1991, L A 9 / J , [14^1, [1JU1 


48 

to 


v^U 


S 41 1 600 


2.4-165 


173 


TS81 T6^1 T6S1 T7S1 T761 T8^1 T871 TQ11 TQ41 ni^l T11S1 T1911 ri9^1 T1941 T1971 
[Joj, L^^j — [UJj, [/jj, [/uj, L^9j — [o/j, ["ij, ["H-j, [iijj, [iijj, [izij, [izj)j, [izh-j, [iz/j, 

ri9Qi n^n n^^i n^si n^7i ri48i nsoi 

[IZ^l, [1911— [1331, L 1 931— [19/1, [14ol— [13U1 


49 


In 


S 41 1 ^09 

J .4-1 — 1 J)UZ 


1 09 16S 
1 .uz— lOJ 


6^ 


T7S1 T8S1 T871 T10S1 T1 1 SI T1971 n^1l Tl^l TI^SI n^71 T1 ^Q1 T1481 T14Q1 

[/JJ, [OJJ — [O/J, [1UJJ, [11JJ, [1Z/J, [131J, [IJJJ, [1 JJJ — [1 J 1 J, [IJ^J, [140J, [147j 


50 


Sn 


S 41 1 408 1 

J.H-l — 14UO. 1 


1.03-165 


405 


TSQ1 T611 T641 T6S1 T711 T761 T811 T8S1 T871 TQ91 TQ^I TQ61 TQ81 TQQ1 T10Q1 
[jyj— [Oij, [041, [OJi, [/ij— L/oj, L°9j— Lo/j, iyz\, L"9i, L"oi, L"ol, L""J, L 1 ^"!, 

T1191 T1161 T11Q1 T1971 n^11 Tl^l TI^SI n^71 T1411 T14S1 T1471 T1 SOI T1 S91 

[11ZJ— [llOJ, [1 l!7j— [1Z / J, [131J, [1JJ>J, [1 J>JJ— [1 D 1 J, [141J, [14JJ, [14/J— [1JUJ, [1JZJ, 

L1S^1 L1SS1 


51 


Sb 


c 41 en <4 
9.41—9:7.94 


QO 1^1 
y\J— 19 1 


94 
Z4 


r7Si n^n n^^i n^6i n^7i ri4Qi 
L/9J, Yl~>l\, Yl~>~>\, LijoJ, lL3/\, iL^y] 


S9 

JZ 


1C 


22.1-59.54 


121-133 


2 


L1111 L14Q1 

[111J, [14^J 


53 


j 


661 6 1 1 90 S 
OO 1 .0— 1 1ZU.J 


QS 110 

7 J— 1 1U 


4 


TSQI T81I 


54 


Xe 


SQ S4 661 6 

J". Jt — UU 1 .u 


^09 190 

J.UZ — 1ZU 


25 


TQS1 T1S11 
L"9J, [1911 


55 


Cs 


17.44 


90 


I 


L1 ^1 1 
[131J 


S6 


r>d 


1 7 44 1 1 90 S 
1 / ,tt — 1 1ZU.J 


57.5-133 


15 


T6^1 T7S1 T811 T1111 T1191 n^11 TL^SI T14Q1 
[091, [/9J, [oil, [iiil, L a A^J , LA911, L 1 9Jl, L a^-^J 


57 


La 


M AA ^ 
1 / .44— OOl .0 


1 77 QO 
1. / l—y\J 


9 


n^ii n^si n/i7i 
L191J, LI9jJ, Ll^'J 


58 


Ce 


14.93-59.54 


90-131 


6 


[75], [131], [132], [135] 


59 


Pr 


14 cn 97Q 9 

lt."J — Z / :7.Z 


90-90 


10 


n9^i n94i n^ii n^^i n^si nsoi 

[IZJJ, [1ZH-J, [131J — [133J, [13JJ, [1JL/J 


60 


Nd 


22.16-661.6 


1.77-135 


11 


[75], [147], [152] 


62 


Sm 


14.93-279.2 


90-133 


13 


[75]! [HI]! [115], [123], [124], [132], [133], [135], [150] 


63 


Eu 


17.44 


90 


1 


[131] 


64 


Gd 


17.44-320 


10-160 


34 


[62], [75], [112], [113], [115], [123], [124], [130], [131], [135], [148]-[150] 


65 


Er 


59.54-320 


45-121 


11 


[130], [149] 


66 


Dy 


16.58-661.6 


1.77-160 


36 


[62], [75], [113], [115], [123], [124], [127], [130], [131], [133], [135], [147]-[150] 


67 


Ho 


16.58-279.2 


10-160 


20 


[75], [111]-[113], [123], [124], [131], [133], [135], [149], [150] 


68 


Er 


16.58-59.54 


10-165 


22 


[62], [113], [115], [121], [131], [133], [135], [148], [149] 


69 


Tm 


22.1 


133 


1 


[HI] 


70 


Yb 


16.58-320 


10-160 


35 


[62], [75], [111]-[113], [115], [123], [124], [130], [131], [133], [135], [149], [150], 
[156] 


72 


Hf 


22.16-511 


0.5-160 


13 


[62], [75], [102], [113] 
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TABLE IV 



Summary of the experimental data used in the validation analysis: atomic numbers 73-92 



Element 


Energy range 


Angle range 


Sample 


References 


Z 


Symbol 


(Ke v ; 


(degrees) 


size 




/ J 


Id 


13.95-1408.1 


0.5-165 


295 


[62], [66], [67], [69], [75], [77], [85]-[87], [98], [102], [105], [112]-[114], [119], [120], 
no^i N97i noon nzKi nzLsi n^m n^9i n^^i 

[1ZJJ— [1Z/J, [IZ^J, [l^JJ, Ll"4-oJ— [1 JUJ, [1JZJ, [1JJJ, [1JOJ 


1A 


W 
vv 


n QC 1 Af)Q 1 


U.J- lOJ 


i ^ i 

i j i 


r^Qi r^m r78i rs^i rs7i nrm ni9i ni^i no^i nozn \mi~\ nasi 
[J^j, L^uj, [OJJ— [OJJ, L'OJ? L"3J — L"'J? L^zj, L^^^J? L^^^J? Li^J? L^^^J? L^^'J? L^^oJ - 

[IjUJ, 


75 


Re 


22.16-511 


0.5-121 


3 


[102], [148], [149] 


77 

/ / 


Ir 


22.16-511 


0.5-160 


11 


[102], [113], [148], [149] 


78 


Pt 


J.^l-l 1ZU.J 


U.J- lOU 


i nn 

1UU 


rozn nrm ni9i n i^i niQi n9m ri97i n^n n^^i n^si n^7i n^si 
[OJj, L-' u> L^uzj, [iizj, [iijj, [iiyj, [izuj, l ± / j , [i JiJ- JjJ? L^^^J? L^-^'J? L^-^"J? 

ri/iO"| MzLQI 


79 


Au 


j.^i-j 1 1 


u.j-iou 


1 08 
1U5 


r^^i r^Qi r^m tqqi yqii nnoi rinii nioi nni ni^i n^ii n^^i 

Looj, ioy\, L/uJ? L'JJ? Lo-^J? liuZ], [lu/j, [iizj, [iijj, [iijj, [ijij— [ijjj, [U/j— 

n<ifn nzisi nzifvi nzisi nziQi ns^i 

[1^+UJ, [l^tJJ, [I'H-OJ, [l^toj, [l^t^J, [1JOJ 


80 


Hg 


22.1-1600 


0.5-135 


16 


[58], [59], [68], [81], [102], [111], [112], [144], [147], [149] 


81 


Tl 


22.1-511 


0.5-133 


3 


[102], [111], [112] 


82 


Pb 


5.41-1613 


0.5-160 


797 


[3], [59]-[61], [64]-[79], [81]-[84], [89], [91]-[94], [96]-[102], [104], [107]-[110], 
[112]-[117], [119], [120], [122]-[129], [131]-[135], [137]-[140], [142]-[150], [152], 
[153], [155], [156], [158] 


83 


Bi 


22.16-511 


0.5-131 


9 


[70], [75], [102], [107], [112], [148], [149] 


90 


Th 


22.16-661.6 


0.5-125 


11 


[102], [112], [123], [124], [147]-[150] 


92 


U 


22.16-1613 


0.5-160 


262 


[63], [75], [92], [93], [100]-[102], [106], [110], [112]-[114], [118], [123], [124], [148]- 



[150], [155] 



digitized cross section values is hindered by the difficulty 
of appraising the experimental points and their error bars in 
plots that may span several orders of magnitude in logarithmic 
scale; these digitized data were not used in the validation 
analysis, if found incompatible with other measurements in 
the same experimental configuration (target element, incident 
photon energy and measured angle) reported in the literature 
in numerical form. 

Large discrepancies are evident in some of the experimental 
data; systematic effects are likely present in some cases, where 
the Wald-Wolfowitz test [57] detects sequences of positive or 
negative differences between data samples originating from 
different experimental groups, which are incompatible with 
randomness. Experimental data exhibiting large discrepancies 
with respect to other measurements in similar configurations, 
that hint to the presence of systematic effects, are excluded 
from the validation process. 

Correct estimate of experimental errors is a concern in the 
validation process, since unrealistic estimation of the experi- 
mental errors may lead to incorrect conclusions regarding the 
rejection of the null hypothesis in tests whose statistic takes 
into account experimental uncertainties explicitly. Although 
technological developments have contributed to improved pre- 
cision of measurement, it has been remarked that in some 
circumstances estimates of experimental uncertainties of the 
order of 5% or less reported in the literature may be optimistic 
[5]. Experimental measurements claiming much smaller uncer- 
tainties than similar ones have been critically evaluated in the 
analysis process; effects related to the validation of simulation 
models are discussed in section IV-BI 

D. Data analysis method 

The evaluation of the simulation models performed in this 
study has two objectives: to validate them quantitatively, and 
to compare their relative capabilities. 

The scope of the software validation process is defined 
according to the guidelines of the pertinent IEEE Standard 



[161], which conforms to the life cycle process standard 
defined in ISO/IEC Standard 12207 [162]. For the problem 
domain considered in this paper, the validation process pro- 
vides evidence that the software models photon-atom elastic 
scattering consistent with experiment. 

A quantitative analysis, based on statistical methods, is 
practically possible only for the validation of differential cross 
sections, for which a large sample of experimental data is 
available. The analysis of differential cross sections is articu- 
lated over two stages: the first one estimates the compatibility 
between the values calculated by each simulation model and 
experimental data, while the second exploits the results of the 
first stage to determine whether the various models exhibit any 
significant differences in their compatibility with experiment. 

The first stage encompasses a number of test cases, each 
one corresponding to a photon energy for which experimental 
data are available. For each test case, cross sections calculated 
by the software are compared with measured ones by means of 
goodness-of-fit tests; the null hypothesis in the tests is defined 
as the equivalence of the simulated and experimental data dis- 
tributions subject to comparison. The Statistical Toolkit [163], 
[164] is exploited for these tests. The level of significance of 
the tests described in the following sections is 0.01, unless 
differently stated; it means that the null hypothesis is rejected 
whenever the p-value resulting from the test statistic is smaller 
than 0.01. 

The goodness-of-fit analysis is primarily based on the x 2 
test [165]. Among goodness-of-fit tests, this test has the 
peculiarity of taking into account experimental uncertainties 
explicitly; therefore the test statistic is sensitive to their correct 
appraisal. Three tests based on the empirical distribution 
function, the Kolmogorov-Smirnov [166], [167], Anderson- 
Darling [168], [169] and Cramer-von Mises [170], [171] tests, 
have also been applied; nevertheless, as documented in the 
following sections, they exhibit limited discriminant power 
over the data sample of this study. In this context one may 
want to recall that the power of goodness-of-fit tests is still 
object of active research in statistics. 
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The "efficiency" of a physics model is defined as the 
fraction of test cases in which the x 2 test does not reject 
the null hypothesis at 0.01 level of significance: it quantifies 
the capability of that simulation model to produce results 
statistically consistent with experiment over the whole set of 
test cases, which in physical terms means over the whole 
energy range of the experimental data sample involved in the 
validation process. 

The second stage of the statistical analysis quantifies the 
differences of the simulation models in compatibility with 
experiment. It consists of a categorical analysis based on 
contingency tables, which derive from the results of the x 2 
test: the outcome of this test is classified as "fail" or "pass", ac- 
cording respectively to whether the hypothesis of compatibility 
of experimental and calculated data is rejected or not. The 
categorical analysis takes as a reference the simulation model 
exhibiting the largest efficiency at reproducing experimental 
data; the other models are compared to it, to determine whether 
they exhibit statistically significant differences of compatibility 
with measurements. 

The null hypothesis in the analysis of a contingency table 
assumes the equivalent compatibility with experiment of the 
models it compares. Contingency tables are analyzed with 
Fisher's exact test [172] and with the x 2 test applying Yates' 
continuity correction [173]; the latter ensures meaningful 
results when one or more cells of the table have few entries. 
Pearson's x 2 test [174] is also performed on contingency 
tables, when their content is consistent with its applicability. 
The use of different tests in the analysis of contingency tables 
contributes to the robustness of the results, as it mitigates the 
risk of introducing systematic effects, which could be due to 
the peculiar mathematical features of a single test. 

The significance level for the rejection of the null hypothesis 
in the analysis of contingency tables is 0.05, unless differently 
specified. 

The scarcity of experimental measurements of elastic scat- 
tering total cross sections in the literature hinders the valida- 
tion of simulation models through similar statistical analysis 
methods: only qualitative general considerations can be made, 
while quantitative assessments are limited by the small sample 
of published experimental data. Nevertheless, the results of 
the analysis of differential cross sections contribute indirectly 
to assess the validity of total cross sections, since some 
of the implemented total cross section models derive from 
the integration of the differential cross sections subject to 
validation. 

V. Results 

The following sections report the results of the comparison 
of differential and total cross sections with experimental data, 
and of the evaluation of the computational performance of the 
various simulation models. 

A. General features 

Elastic scattering cross sections vary by several orders of 
magnitude over the data sample involved in the validation anal- 
ysis. Cross sections decrease as a function of photon energy 



and scattering angle, and are larger for heavy elements: these 
trends are illustrated in Figs. [T] and [2] that show experimental 
photon elastic scattering measurements along with the values 
calculated by representative simulation models for the same 
energy and scattering angle settings. 

Experimental differential cross sections are the result of 
all the physics processes that contribute to photon elastic 
scattering, while the simulation models evaluated in this paper 
account for the Rayleigh scattering amplitude only or, in the 
case of the model based on S -matrix calculations, for the 
sum of Rayleigh and Thomson scattering amplitudes. This 
feature is evident in Fig. [3] which includes some of the higher 
energy measurements in the experimental data sample: other 
processes, such as Delbruck scattering, should be taken into 
account in the simulation, along with Rayleigh scattering, to 
model photon elastic scattering accurately at higher energies. 
The plots also expose some characteristics of the experimental 
data: systematic effects affecting some of the measurements, 
and the presence of outliers in the experimental sample. 

The distributions in Figs. [4] and [5] illustrate the difference 
between calculated and measured differential cross sections, 
for a few representative models: two options based on the 
form factor approximation, respectively using the form fac- 
tors tabulated in EPDL97 and modified form factors with 
anomalous scattering factors, and the model based on S -matrix 
calculations. Fig. [4] shows the relative difference between 
simulated and experimental values. Fig. [5] accounts for the 
different precision of the data in the experimental sample: it 
represents the difference between calculated and experimental 
differential cross sections in units of standard deviations, i.e. 
of the uncertainties associated with the experimental mea- 
surements. These plots are limited to photon energies up to 
1.5 MeV, since, as shown in Fig. [3] at higher energies other 
processes should be taken into account to describe photon 
elastic scattering. The long tails of the distributions in Figs. [4] 
and [5] are the result of multiple factors: neglect of amplitudes 
other than Rayleigh scattering in some of the calculations, 
inadequacy of the simulation models, and contamination of the 
experimental sample by measurements affected by systematic 
effects and outliers. One can observe in Fig. [4] that the 
distribution of the relative differences between the differential 
cross sections calculated by the models and experimental data 
is narrower for the S-matrix model, while it is wider for the 
cross sections based on EPDL and on modified form factors 
with anomalous scattering factor corrections. 

The one sample t-test [175] for the hypothesis of null mean 
difference between simulation and experiment confirms the 
qualitative features observed in Figs. [4] and [5] Fig. [6] shows 
the 99% confidence intervals for the true mean difference 
associated with the various models: the model based on S- 
matrix calculations is the only one for which the confidence 
interval includes zero. The hypothesis of null mean difference 
from experiment equal to zero is rejected for all the cross 
section models with 0.001 significance, with the exception of 
the model based on S-matrix calculations, for which the p- 
value resulting from the t-test is 0.500. 

The algorithm for Rayleigh scattering simulation reengi- 
neered from Penelope 2008-2011 versions produces very 
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Fig. 1. Differential cross section as a function of energy at 9 = 90° for aluminium, copper, tin and lead: experimental measurements (black circles), 
calculations based on S-matrix (SM, black empty squares) and on EPDL (red circles). The sources of experimental data are documented in Tables [Till and |IV| 



similar results to the algorithm originally developed in [45], 
which corresponds to the policy class identified in this paper 
as EPDL: in fact, both of them are based on EPDL97 data. 
Figs. [7] and [8] show the relative difference between the total and 
differential cross sections calculated by these two models; the 
small differences are due to the different reference grids used 
by the two models in the respective tabulations of total cross 
sections and non-relativistic form factors, which affect the 
values calculated by interpolation. The difference of resulting 
total cross sections appears insignificant, considering that the 
scale of Fig. [7] extends to relative differences up to 0.05%, 



while the distribution of differences between differential cross 
sections, shown in Fig. [8] is wider, having a standard deviation 
of 1%: the effect of these differences on the accuracy of the 



two models is quantitatively estimated in section V-B 



B. Differential cross sections 

The analysis of differential cross sections is articulated over 
the set of discrete energies for which experimental data are 
available. The p-values resulting from the \ 2 test applied 
to the data samples corresponding to each energy are listed 
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Fig. 2. Differential cross section as a function of scattering angle for representative energies and target elements: experimental measurements (black circles), 
calculations based on S-matrix (SM, black empty squares) and on EPDL (red circles). The S-matrix calculations account for Rayleigh scattering and nuclear 
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TABLE V 

Differential corss sections: p- values resulting from the x 2 test over the whole angular range 



E 


Penelope 


Penelope 


EPDL 


EPDL 


SM 


RF 


NF 


MF 


MF 


RF 


(keV) 


2001 


2008-2011 




ASF 










ASF 


ASF 


5.45 


0.998 


0.996 


0.992 


0.297 


1.000 


0.952 


0.996 


0.999 


1.000 


1.000 


5.93 


0.706 


0.745 


0.741 


0.049 


0.999 


0.432 


0.732 


0.775 


1.000 


1.000 


6.4 


0.860 


0.830 


0.863 


0.148 


1.000 


0.728 


0.829 


0.923 


1.000 


1.000 


6.93 


0.913 


0.877 


0.876 


0.008 


1.000 


0.722 


0.883 


0.944 


1.000 


1.000 


7.47 


0.695 


0.626 


0.614 


0.041 


1.000 


0.386 


0.637 


0.749 


1.000 


1.000 


8.04 


0.053 


0.001 


0.001 


< 0.001 


0.696 


< 0.001 


0.001 


< 0.001 


0.736 


0.727 


8.63 


0.144 


0.249 


0.241 


< 0.001 


0.759 


0.020 


0.199 


0.232 


0.820 


0.812 


9.23 


0.023 


0.033 


0.030 


< 0.001 


0.652 


0.002 


0.024 


0.029 


0.719 


0.708 


13.95 


< 0.001 


< 0.001 


< 0.001 


0.002 


0.154 


< 0.001 


< 0.001 


< 0.001 


0.046 


0.064 


14.93 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


16.58 


0.956 


0.908 


0.901 


0.587 


0.218 


1.000 


0.947 


0.947 


0.305 


0.281 


17.44 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.102 


< 0.001 


< 0.001 


< 0.001 


0.013 


0.025 


17.75 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.105 


< 0.001 


< 0.001 


< 0.001 


0.009 


0.025 


21.12 


0.008 


0.001 


0.001 


0.017 


0.099 


< 0.001 


0.001 


0.001 


0.124 


0.118 


22.1 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.014 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


23.1 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.002 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


24.14 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.002 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


25.2 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


26.27 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.001 


< 0.001 


< 0.001 


0.001 


< 0.001 


26.36 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.772 


0.758 


30.85 


< 0.001 


0.280 


0.322 


< 0.001 


0.011 


0.955 


0.399 


0.507 


0.028 


0.015 


30.97 


0.001 


0.011 


0.020 


0.004 


0.377 


0.001 


0.006 


0.002 


0.667 


0.509 


32.06 


0.003 


0.078 


0.100 


0.197 


0.602 


0.030 


0.066 


0.036 


0.885 


0.749 


33.3 


0.004 


0.077 


0.082 


0.048 


0.306 


0.154 


0.104 


0.072 


0.589 


0.422 


34.6 


< 0.001 


< 0.001 


< 0.001 


0.005 


0.033 


0.006 


0.001 


0.001 


0.071 


0.046 


35.86 


0.001 


0.001 


0.001 


0.008 


0.148 


0.195 


0.002 


0.002 


0.314 


0.179 


39.9 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.061 


0.043 


0.004 


0.003 


0.083 


0.072 


42.75 


< 0.001 


0.001 


0.001 


0.002 


0.183 


0.001 


0.003 


0.002 


0.082 


0.089 


47.24 


< 0.001 


0.011 


0.021 


0.003 


0.019 


0.005 


0.002 


0.065 


0.022 


0.044 


59.54 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


81 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


84.3 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


88 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


122 


0.196 


0.500 


0.541 


< 0.001 


0.107 


< 0.001 


0.478 


0.287 


0.033 


0.182 


136 


0.534 


0.837 


0.790 


< 0.001 


0.913 


0.047 


0.893 


0.836 


0.669 


0.891 


145.4 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.037 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


244.7 


< 0.001 


0.025 


0.032 


0.006 


0.038 


< 0.001 


0.027 


0.024 


0.037 


0.036 


279.2 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.201 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


317 


< 0.001 


< 0.001 


< 0.001 


0.019 


< 0.001 


0.200 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


320 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


334 


0.031 


0.873 


0.801 


< 0.001 


0.118 


< 0.001 


0.617 


0.757 


< 0.001 


0.022 


344 


< 0.001 


0.047 


0.033 


< 0.001 


0.080 


< 0.001 


0.040 


0.140 


< 0.001 


< 0.001 


411.8 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.136 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


443.96 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.053 


< 0.001 


< 0.001 


0.111 


0.012 


< 0.001 


465 


0.758 


0.542 


0.518 


0.488 


0.285 


0.120 


0.535 


0.302 


0.287 


0.290 


468.1 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.053 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


511 


0.754 


0.752 


0.756 


0.756 


0.730 


0.746 


0.753 


0.727 


0.731 


0.731 


661.6 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.044 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


723 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.096 


< 0.001 


< 0.001 


0.114 


0.009 


< 0.001 


779 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.683 


< 0.001 


< 0.001 


0.333 


< 0.001 


< 0.001 


811 


4.561 


6.060 


6.245 


6.295 


5.822 


8.857 


6.174 


5.828 


5.881 


5.784 


847 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.598 


< 0.001 


< 0.001 


0.870 


0.278 


< 0.001 


868 


0.005 


0.162 


0.141 


0.043 


0.919 


< 0.001 


0.142 


0.925 


0.854 


< 0.001 


878 


0.026 


0.140 


0.119 


0.115 


0.164 


0.001 


0.123 


0.153 


0.140 


0.156 


889 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.323 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


952 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


964 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.152 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1005 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.987 


< 0.001 


< 0.001 


0.968 


0.697 


< 0.001 


1019 


0.280 


0.053 


0.045 


0.045 


0.075 


< 0.001 


0.049 


0.073 


0.066 


0.095 


1038 


< 0.001 


0.041 


0.011 


0.011 


0.004 


< 0.001 


0.056 


0.019 


< 0.001 


< 0.001 


1086 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.125 


< 0.001 


< 0.001 


0.011 


< 0.001 


< 0.001 


1112 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.060 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1120 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1173 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.029 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1189 


0.948 


0.401 


0.380 


0.380 


0.523 


0.022 


0.384 


0.482 


0.468 


0.522 


1238 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.011 


< 0.001 


< 0.001 


0.008 


< 0.001 


< 0.001 


1274 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.076 


< 0.001 


< 0.001 


0.016 


< 0.001 


< 0.001 


1302 


0.054 


0.039 


0.031 


0.031 


0.050 


< 0.001 


0.032 


0.048 


0.044 


0.057 


1332 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1359 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.382 


< 0.001 


< 0.001 


0.325 


0.133 


< 0.001 


1408 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.312 


< 0.001 


< 0.001 


0.068 


< 0.001 


< 0.001 


>1600 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 
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TABLE VI 

Differential cross sections: p- values resulting from the x 2 test, data sample up to 90° 



E 


Penelope 


Penelope 


EPDL 


EPDL 


SM 


RF 


NF 


MF 


MF 


RF 


(keV) 




2001 


2008-2011 




ASF 










ASF 


ASF 


5.45 




0.998 


0.996 


0.992 


0.297 


1.000 


0.952 


0.996 


0.999 


1.000 


1.000 


5.9 




0.706 


0.745 


0.741 


0.049 


0.999 


0.432 


0.732 


0.775 


1.000 


1.000 


6.40 




0.860 


0.830 


0.863 


0.148 


1.000 


0.728 


0.829 


0.923 


1.000 


1.000 


6.93 




0.913 


0.877 


0.876 


0.008 


1.000 


0.722 


0.883 


0.944 


1.000 


1.000 


7.47 




0.695 


0.626 


0.614 


0.041 


1.000 


0.386 


0.637 


0.749 


1.000 


1.000 


8.04 




0.053 


0.001 


0.001 


< 0.001 


0.696 


< 0.001 


0.001 


< 0.001 


0.736 


0.727 


8.63 




0.144 


0.249 


0.241 


< 0.001 


0.759 


0.020 


0.199 


0.232 


0.820 


0.812 


9.23 




0.023 


0.033 


0.030 


< 0.001 


0.652 


0.002 


0.024 


0.029 


0.719 


0.708 


14.93 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


16.58 




0.956 


0.908 


0.901 


0.587 


0.218 


1.000 


0.947 


0.947 


0.305 


0.281 


17.4 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.102 


< 0.001 


< 0.001 


< 0.001 


0.013 


0.025 


21.12 




0.008 


0.001 


0.001 


0.017 


0.099 


< 0.001 


0.001 


0.001 


0.124 


0.118 


22.1 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.159 


< 0.001 


< 0.001 


< 0.001 


0.021 


0.040 


23.10 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.002 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


24.14 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.002 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


25.2 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


26.27 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.001 


< 0.001 


< 0.001 


0.001 


< 0.001 


30.9 


< 


0.001 


0.280 


0.322 


< 0.001 


0.011 


0.955 


0.399 


0.507 


0.028 


0.015 


31.0 




0.001 


0.011 


0.020 


0.004 


0.377 


0.001 


0.006 


0.002 


0.667 


0.509 


32.06 




0.003 


0.078 


0.100 


0.197 


0.602 


0.030 


0.066 


0.036 


0.885 


0.749 


33.30 




0.004 


0.077 


0.082 


0.048 


0.306 


0.154 


0.104 


0.072 


0.589 


0.422 


34.60 


< 


0.001 


< 0.001 


< 0.001 


0.005 


0.033 


0.006 


0.001 


0.001 


0.071 


0.046 


35.86 




0.001 


0.001 


0.001 


0.008 


0.148 


0.195 


0.002 


0.002 


0.314 


0.179 


39.9 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.061 


0.043 


0.004 


0.003 


0.083 


0.072 


42.75 


< 


0.001 


0.001 


0.001 


0.002 


0.183 


0.001 


0.003 


0.002 


0.082 


0.089 


47.24 


< 


0.001 


0.011 


0.021 


0.003 


0.019 


0.005 


0.002 


0.065 


0.022 


0.044 


59.54 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


0.001 


81 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


84.3 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


122 




0.196 


0.500 


0.541 


< 0.001 


0.107 


< 0.001 


0.478 


0.287 


0.033 


0.182 


136 




0.534 


0.837 


0.790 


< 0.001 


0.913 


0.047 


0.893 


0.836 


0.669 


0.891 


145.4 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.459 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


244.7 


< 


0.001 


0.025 


0.032 


0.006 


0.038 


< 0.001 


0.027 


0.024 


0.037 


0.036 


279.2 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.096 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


317 


< 


0.001 


< 0.001 


< 0.001 


0.019 


< 0.001 


0.200 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


320 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


1.000 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


334 




0.031 


0.873 


0.801 


< 0.001 


0.118 


< 0.001 


0.617 


0.757 


< 0.001 


0.022 


344 


< 


0.001 


0.047 


0.033 


< 0.001 


0.080 


< 0.001 


0.040 


0.140 


< 0.001 


< 0.001 


411.8 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.385 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


443.96 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.053 


< 0.001 


< 0.001 


0.111 


0.012 


< 0.001 


465 




0.758 


0.542 


0.518 


0.488 


0.285 


0.120 


0.535 


0.302 


0.287 


0.290 


468.1 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.051 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


511 




0.754 


0.752 


0.756 


0.756 


0.730 


0.746 


0.753 


0.727 


0.731 


0.731 


661.6 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.028 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


723 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.096 


< 0.001 


< 0.001 


0.114 


0.009 


< 0.001 


779 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.683 


< 0.001 


< 0.001 


0.333 


< 0.001 


< 0.001 


811 




0.102 


0.048 


0.044 


0.043 


0.054 


0.012 


0.046 


0.054 


0.053 


0.055 


847 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.598 


< 0.001 


< 0.001 


0.870 


0.278 


< 0.001 


868 




0.005 


0.162 


0.141 


0.043 


0.919 


< 0.001 


0.142 


0.925 


0.854 


< 0.001 


878 




0.026 


0.140 


0.119 


0.115 


0.164 


0.001 


0.123 


0.153 


0.140 


0.156 


889 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.042 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


952 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


964 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.152 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1005 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.987 


< 0.001 


< 0.001 


0.968 


0.697 


< 0.001 


1019 




0.280 


0.053 


0.045 


0.045 


0.075 


< 0.001 


0.049 


0.073 


0.066 


0.095 


1038 


< 


0.001 


0.041 


0.011 


0.011 


0.004 


< 0.001 


0.056 


0.019 


< 0.001 


< 0.001 


1085.87 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.125 


< 0.001 


< 0.001 


0.011 


< 0.001 


< 0.001 


1112 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.060 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1120 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.090 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1173 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.090 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1189 




0.948 


0.401 


0.380 


0.380 


0.523 


0.022 


0.384 


0.482 


0.468 


0.522 


1238 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.011 


< 0.001 


< 0.001 


0.008 


< 0.001 


< 0.001 


1274 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.076 


< 0.001 


< 0.001 


0.016 


< 0.001 


< 0.001 


1302 




0.054 


0.039 


0.031 


0.031 


0.050 


< 0.001 


0.032 


0.048 


0.044 


0.057 


1332 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


1359 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.382 


< 0.001 


< 0.001 


0.325 


0.133 


< 0.001 


1408 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


0.312 


< 0.001 


< 0.001 


0.068 


< 0.001 


< 0.001 


>1600 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 
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Fig. 3. Differential cross section as a function of energy at = 90° for tin, tantalum, lead and uranium: experimental measurements (black circles), 
calculations based on S-matrix (SM, black empty squares) and on EPDL (red circles). The sources of experimental data are documented in Tables [Till and |IV| 



in Tables |VJ |Vl] and VII respectively concerning the whole 
angular range, data up to 90° and above 90° (backscattering). 

The efficiency of the various simulation models is reported 
in Table |VIII| The results are listed for three data samples: 
the whole experimental data collection, the data for scattering 
angles up to 90° and above 90°. The simulation model based 
on the interpolation of RTAB S-matrix calculations, indicated 
in the tables as SM, stands out as the one exhibiting the largest 
efficiency in all test configurations. 

Three test cases for which the x 2 test rejects the hypothesis 
of compatibility between the S-matrix model and experimental 



data involve measurements from a single reference; due to 
the lack of other experimental data in the same configuration, 
one cannot exclude that the incompatibility between S-matrix 
calculations and experiment could be ascribed to systematic 
effects in the data, rather than to deficiencies of the simulation 
model. The exclusion of these three test cases from the 
analysis slightly increases the efficiencies of all the simulation 
models (e.g. 0.81±0.05 for the S-matrix model and 0.40±0.05 
for the EPDL-based model, over the whole angular range), 
nevertheless it does not substantially change the conclusions 
one can draw from the x 2 test. 
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TABLE VII 

Differential cross sections: p- values resulting from the x 2 test, data sample above 90° 



E 

(keV) 


Penelope 
2001 


Penelope 
2008-2011 


EPDL 


EPDL 

ASF 


SM 




RF 




NF 






MF 




MF 
ASF 




RF 

ASF 


13.95 


< 


n nni 
U.UU1 


< 


n nni 
U.UUl 


< 0.001 


0.002 


0.154 


< 


0.001 


< 


0.001 


< 


o 


001 




0.046 




0.064 


17.75 


< 


n nni 
U.UUl 


< 


n nni 
U.UUl 


< 0.001 


< 0.001 


0.105 


< 


0.001 


< 


0.001 


< 


o 


001 




0.009 




0.025 


22.1 


< 


n nni 
U.UUl 


< 


n nni 
U.UUl 


< 0.001 


< 0.001 


0.010 




0.001 




0.001 




o 


001 




0.001 




0.001 


26.36 


< 


n nni 

U.UUl 


< 


n nni 

U.UUl 


< 0.001 


< 0.001 


< 0.001 


< 


0.001 


< 


0.001 


< 


o 


001 




0.772 




0.758 


59.54 


< 


n nni 
U.UUl 


< 


n nni 
U.UUl 


< 0.001 


< 0.001 


< 0.001 


< 


0.001 


< 


0.001 


< 


o 


001 


< 


0.001 




0.001 


81 




0.001 




0.001 


< 0.001 


< 0.001 


0.365 


< 


0.001 


< 


0.001 


< 


o 


001 




0.004 




0.082 


88 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 


0.001 


< 


0.001 


< 


o 


001 


< 


0.001 


< 


0.001 


145.4 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


279.2 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


0.096 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


317 


< 


0.001 


< 


0.001 


< 0.001 


0.019 


< 0.001 




0.200 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


411.8 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


0.034 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


468.1 




0.015 




0.047 


0.037 


< 0.001 


0.368 


< 


0.001 




0.045 


< 





001 


< 


0.001 


< 


0.001 


661.6 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


0.508 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


889 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


0.940 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


1120 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


1173 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


0.022 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 


1332 


< 


0.001 


< 


0.001 


< 0.001 


< 0.001 


< 0.001 


< 


0.001 


< 


0.001 


< 





001 


< 


0.001 


< 


0.001 




Fig. 4. Relative difference between the differential cross sections calculated 
by simulation models and the corresponding experimental measurements: 
calculations based on based S-matrix (SM, solid black line), EPDL (dashed 
red line) and MFASF form factors (dot-dashed blue line). The plot is limited 
to photon energies up to 1.5 MeV 



Fig. 5. Difference between the differential cross sections calculated by sim- 
ulation models and the corresponding experimental measurements, expressed 
in terms of number of standard deviations (experimental uncertainties): 
calculations based on based S-matrix (SM, solid black line), EPDL (dashed 
red line) and MFASF form factors (dot-dashed blue line). The plot is limited 
to photon energies up to 1.5 MeV 



Several cases of inconsistency between simulation models 
and experimental data concern measurements at photon en- 
ergies close to the K or L shell electron binding energies 
of target atom; a number of them involve measurements for 
which relatively small experimental uncertainties are reported, 
that could have been underestimated. The same analysis has 
also been performed on a sample that excludes these data; 
its results, which are reported in T able | IX] confirm the same 
general trend observed in Table |VIIl[ The larger relative 
increase in efficiency of the models based on RTAB observed 
in this analysis with respect to the effect on other models 
suggests that this database might have adopted a coarse grid 
for its tabulations, that degrades the accuracy of interpolation 
in sensitive areas. 



The compatibility between simulation models and experi- 
mental data has been evaluated also through the Kolmogorov- 
Smirnov, Anderson-Darling and Cramer-von Mises goodness- 
of-fit tests. The results, listed in Table |X| suggest that these 
tests exhibit scarce sensitivity to differences between the simu- 
lation models and the experimental data subject to comparison. 

The relative accuracy of the various simulation models is 
compared by contingency tables based on the outcome of the 
X 2 test. In these comparisons the model based on the S-matrix, 
which exhibits the largest efficiency, is taken as a reference; 
contingency tables examine whether the results of the x 2 test 
over other models are statistically equivalent to those of this 
model. 

Since several general purpose Monte Carlo systems base 
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TABLE VIII 

X 2 TEST OUTCOME: TEST CASES COMPATIBLE WITH EXPERIMENT AT 0.01 SIGNIFICANCE LEVEL 



Test 




Penelope 


Penelope 


¥7 niAi 

iLrUL, 


iLrDL, 




Kr 


TVTC 

JMr 


Mr 


A AT? 

Mr 1 


Kr 






2001 


2008-2011 




ASF 










ASF 


ASF 




Test cases 


71 


71 


71 


71 


71 


71 


71 


71 


71 


71 




Pass 


19 


27 


27 


18 


55 


18 


25 


35 


37 


34 


all 


Fail 


52 


44 


44 


53 


16 


53 


46 


36 


34 


37 




Efficiency 


U.Z / 




U.jo 


U.Zj 


U. / / 


U.Zj 


U.3j 


n aq 
UAy 


U.jz 


d AQ 

U.4o 




Error 


±0.05 


±0.06 


±0.06 


±0.05 


±0.06 


±0.05 


±0.06 


±0.06 


±0.06 


±0.06 




Test cases 


67 


67 


67 


67 


67 


67 


67 


67 


67 


67 




Pass 


19 


27 


27 


18 


55 


18 


25 


35 


36 


32 


e < 90° 


Fail 


48 


40 


40 


49 


12 


49 


42 


32 


31 


35 




Efficiency 


0.28 


0.40 


0.40 


0.27 


0.82 


0.27 


0.37 


0.52 


0.54 


0.48 




Error 


±0.05 


±0.06 


±0.06 


±0.05 


±0.05 


±0.05 


±0.06 


±0.06 


±0.06 


±0.06 




Test cases 


17 


17 


17 


17 


17 


17 


17 


17 


17 


17 




Pass 


1 


1 


1 


1 


10 


1 


1 





2 


4 


> 90° 


Fail 


16 


16 


16 


16 


7 


16 


16 


17 


15 


13 




Efficiency 


0.06 


0.06 


0.06 


0.06 


0.59 


0.06 


0.06 


< 0.06 


0.12 


0.24 




Error 


±0.06 


±0.06 


±0.06 


±0.06 


±0.12 


±0.06 


±0.06 




±0.08 


±0.10 



TABLE IX 

X 2 TEST OUTCOME, EXCLUDING DATA SENSITIVE TO K AND L SHELL EFFECTS: TEST CASES COMPATIBLE WITH EXPERIMENT AT 0.01 SIGNIFICANCE 

LEVEL 



Test 




Penelope 


Penelope 


EPDL 


EPDL 


SM 


RF 


NF 


MF 


MF 


RF 






2001 


2008-2011 




ASF 










ASF 


ASF 




Test cases 


71 


71 


71 


71 


71 


71 


71 


71 


71 


71 




Pass 


19 


28 


28 


18 


62 


21 


26 


36 


43 


40 


all 


Fail 


52 


43 


43 


53 


9 


50 


45 


35 


28 


31 




Efficiency 


0.27 


0.39 


0.39 


0.35 


0.87 


0.30 


0.37 


0.51 


0.61 


0.56 




Error 


±0.05 


±0.06 


±0.06 


±0.06 


±0.04 


±0.05 


±0.06 


±0.06 


±0.06 


±0.06 




Test cases 


67 


67 


67 


67 


67 


67 


67 


67 


67 


67 




Pass 


19 


28 


28 


25 


61 


21 


26 


36 


42 


38 


< 90° 


Fail 


48 


39 


39 


42 


6 


46 


41 


31 


25 


29 




Efficiency 


0.28 


0.42 


0.42 


0.37 


0.91 


0.31 


0.39 


0.54 


0.63 


0.57 




Error 


±0.06 


±0.06 


±0.06 


±0.06 


±0.04 


±0.06 


±0.06 


±0.06 


±0.06 


±0.06 




Test cases 


17 


17 


17 


17 


17 


17 


17 


17 


17 


17 




Pass 


1 


1 


1 


1 


11 


1 


1 





3 


4 


> 90° 


Fail 


16 


16 


16 


16 


6 


16 


16 


17 


14 


13 




Efficiency 


0.06 


0.06 


0.06 


0.06 


0.65 


0.06 


0.06 


< 0.06 


0.18 


0.24 




Error 


±0.06 


±0.06 


±0.06 


±0.06 


±0.12 


±0.06 


±0.06 




±0.09 


±0.10 



TABLE X 

GOODNESS-OF-FIT TESTS: TEST CASES COMPATIBLE WITH EXPERIMENT AT 0.05 SIGNIFICANCE LEVEL 



Test 




Penelope 


Penelope 


EPDL 


EPDL 


SM 


RF 


NF 


MF 


MF 


RF 






2001 


2008-2011 




ASF 










ASF 


ASF 




Test cases 


71 


71 


71 


71 


71 


71 


71 


71 


71 


71 




Pass 


67 


68 


68 


61 


69 


62 


67 


62 


69 


64 


Kolmogorov 


Fail 


4 


3 


3 


10 


2 


9 


4 


9 


2 


7 


Smirnov 


Efficiency 


0.94 


0.96 


0.96 


0.86 


0.97 


0.87 


0.94 


0.87 


0.97 


0.90 




Error 


±0.03 


±0.02 


±0.02 


±0.04 


±0.02 


±0.04 


±0.03 


±0.04 


±0.02 


±0.04 




Pass 


68 


69 


69 


62 


70 


63 


69 


65 


70 


63 


Anderson 


Fail 


3 


2 


2 


9 


1 


8 


2 


6 


1 


8 


Darling 


Efficiency 


0.96 


0.97 


0.97 


0.87 


0.99 


0.89 


0.97 


0.92 


0.99 


0.89 




Error 


±0.02 


±0.02 


±0.02 


±0.04 


±0.01 


±0.04 


±0.02 


±0.03 


±0.01 


±0.04 




Pass 


69 


69 


69 


62 


70 


66 


69 


67 


70 


63 


Cramer 


Fail 


2 


2 


2 


9 


1 


5 


2 


4 


1 


8 


von Mises 


Efficiency 


0.97 


0.97 


0.97 


0.87 


0.99 


0.93 


0.97 


0.94 


0.99 


0.89 




Error 


±0.02 


±0.02 


±0.02 


±0.04 


±0.01 


±0.03 


±0.02 


±0.03 


±0.01 


±0.04 
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Fig. 6. Confidence intervals for the true mean difference between differential 
cross sections calculated by Rayleigh scattering simulation models and 
experimental measurements; the dashed red line identifies the value of zero 
in the dot chart. 
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Fig. 8. Relative difference between the differential cross sections calculated 
based on interpolation of Penelope 2008-2011 tabulated form factors derived 
from EPDL97 and by interpolation of native EPDL97 tabulations, as described 
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their simulation of photon elastic scattering on EPDL, it is 
interesting to quantify whether the differential cross sections 
based on this data library differ significantly in compatibility 
with experiment from those derived from S-matrix calcula- 



tions. This analysis is summarized in Table XI which reports 



results for the whole data sample, and in Table XII which 
excludes data at energies close to K and L shell binding 
energies; both tables list results calculated over the whole data 
sample and for scattering angles up to, and above 90°. The 
hypothesis of statistically equivalent capability to reproduce 
experimental measurements by the S-matrix simulation model 
and by the model based on EPDL is rejected at 0.01 level of 
significance. 



Tables [XI] and |XII| also report a comparison between the S- 
matrix model and the model exploiting modified form factors 
with anomalous scattering factors (MFASF), which exhibits 



the highest efficiencies in Tables VIII and IX among those 
based on the form factor approximation. The hypothesis of 
equivalent accuracy is rejected with 0.01 significance also 
for the MFASF model over the whole data sample and the 
subset of data with scattering angles up to 90°, while it is not 
rejected for the backscattering data sample, although with p- 
values close to the critical region. A fortiori, simulation models 
exhibiting lower efficiencies than the MFASF model in Table 



XI] are significantly less accurate than the model based on 
RTAB S-matrix calculations. 

From this analysis one can draw the conclusion that the 
simulation model based on RTAB tabulations of S-matrix 
calculations ensures significantly more accurate results than 
all the other simulation alternatives considered in this study, 



18 



including the EPDL data library currently used by several 
general purpose Monte Carlo systems. The model exploiting 
modified form factors with anomalous scattering factors is 
the most accurate among those based on the form factor 
approximation. 

The model based on pristine EPDL97 tabulations of non- 
relativistic form factors and the Penelope 2008-2011 model, 
based on modified tabulations derived from EPDL97, produce 
equivalent differential cross sections: their efficiencies listed in 
Table |Xl| and |XII| are identical, despite the small differences 



TABLE XI 

Contingency tables related to the whole data sample 



visible in Fig. [8 The tabulations of non-relativistic form 
factors included in RTAB produce slightly less accurate cross 
sections than the models based on EPDL97, although the 
differences in efficiencies are compatible within one standard 
deviation. 

The simulation model based on relativistic form factors 
is less accurate than models exploiting non-relativistic form 
factors: it has been previously noted [109], [176] that the use 
of relativistic wavefunctions in the calculation of form factors 
often produces less accurate results than use of nonrelativistic 
wavefunctions, although - to the best of our knowledge - 
the relative efficiency of these two calcluation methods at 
reproducing experimental data has not been yet quantified with 
statistical methods. 

The Penelope 2001 model is less accurate than more recent 
versions of the code. 

The inclusion of anomalous scattering factors in the cal- 
culations based on EPDL97 does not contribute to improve 
compatibility with experiment, while accounting for anoma- 
lous scattering improves the compatibility with experiment of 
calculations exploiting relativistic and modified form factors. 

The simplified model in G4XrayRayleighModel, which im- 
plements elastic scattering from a point-like charge, produces a 
largely unrealistic description of photon elastic scattering: this 
is visible, for instance, in Fig. [9] which shows the distribution 
of the scattering angle of 10 keV photons interacting with 
carbon generated by this model, compared to the distribution 
based on S -matrix calculations. Given its evident inadequacy, 
this model was not included in the statistical analysis for the 
validation of differential cross sections. 

C. Total Cross Sections 

Total cross sections calculated by a set of simulation models 



are shown in Fig. 10 for four representative elements. Some 
of them are associated with modeling approaches whose capa- 
bilities have been evaluated in the previous section regarding 
differential cross sections: the models based on S-matrix 
calculations, EPDL97, modified form factors with anomalous 
scattering factors and Penelope 2001 parameterizations. Others 
(XCOM, Storm and Israel, Brennan and Cowan), are specific 
for total cross section calculation only. In addition, the most re- 
cent development in the field, Geant4 G4XrayRayleighModel, 
is illustrated in the plots. 

The total cross sections calculated by XCOM, Penelope 
2001, Storm and Israel, and G4XrayRayleighModel appear 
insensitive to the underlying atomic structure; for photon 
energies below 10 keV these models overestimate the total 
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Fig. 9. Probability distribution of the scattering angle of 10 keV photons 
interacting with carbon generated by G4XrayRayleighModel (magenta dashed 
line) this model, compared to the distribution based on S-matrix calculations 
(black solid line; both distributions are normalized to 1. The distribution of 
G4XrayRayleighModel corresponds to scattering from a point-like charge. 
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TABLE XII 

Contingency tables excluding energies close to K and L shell 

BINDING ENERGIES 
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cross section with respect to the S-matrix model, which is 
the result of the integration of the differential cross sections 
assessed in section [y-B] as best reproducing experimental data. 
The parameterized cross sections by Brennan and Cowan 
attempt to account for atomic shell effects at low energies, 
but their results appear largely approximate with respect to 
S-matrix calculations. 

Although in the higher energy end all the calculated cross 



sections look quite similar on the scale of Fig. 10 differences 
are visible in Fig. [TT] that highlights the behaviour of total 
cross section models in a region where experimental data are 
available in the literature [159]. 



Most of the models shown in Fig. [TT] produce qualitatively 
similar results in that energy range, apart from the interpolation 
of Storm and Israel's cross sections, the parameterizations by 
Brennan and Cowan, and G4XrayRayleighModel\ this model 
exhibits large discrepancies with respect to Gowda's measure- 
ments. The total elastic scattering cross sections calculated 
by the S-matrix, EPDL and MFASF models are compatible 
with the measurements in [159] with 0.01 significance; the 
cross sections based on Storm and Israel's tabulations and 
calculated by XCOM lie respectively 2.1 and 1.1 standard 



deviations away from the experimental values for barium, 
while they differ from the experimental values of the other 
three target elements by 3.4 to 5.2 standard deviations. The 
cross sections calculated by G4XrayRayleighModel lie 26 to 
41 standard deviations away from the experimental values. 

Experimental data from [160] and total cross sections cal- 
culated by various simulation models in the energy range 



between 250 and 700 keV are in Fig. 12 for 9 elements; the 
experimental errors of the data at 279.2 keV reported in the 
figure have been scaled by a factor 10 with respect to the 
values reported in [160], due to the apparent inconsistency 



discussed in section |IV-C| The plots for the other 9 target 
elements measured in [160] exhibit similar behaviour to those 



shown in Fig. 12 The source of the data [160], which has 
not been subject to peer review for publication in a scientific 
journal, and the concerns about the experimental uncertainties 



discussed in section IV-C suggest caution in using these 
data for quantitative validation of the simulation models; 
nevertheless, these plots, as well as those for the other 9 
measured elements, confirm qualitatively the same conclusions 
drawn from the analysis of the subset of published data. 

Total cross sections based on EPDL97 appear quite similar 
to those deriving from the integration of S-matrix differen- 
tial cross sections, whose accuracy has been quantitatively 



assessed in section V-B In the energy range between 5 keV 
and 1.5 MeV, which approximately corresponds to the scope 
of the validation of differential cross sections documented 
in section |V-B| the Kolmogorov-Smirnov test finds EPDL97 
total cross sections compatible with those based on S-matrix 
calculations for 72% ± 5% of the elements, while the total 
cross sections based on modified form factors with anomalous 
scattering factors are equivalent to those based on S-matrix 
for 97% ± 2% of the test elements. The significance of these 
tests is 0.01. 

XCOM cross sections are found compatible with those 
based on S-matrix in 28% ± 5% of the test cases; this result 
can be considered representative also for a set of models 
(Storm and Israel's, Brennan and Cowan's, Penelope 2001) 
that exhibit similar characteristics. 

Based on these considerations, one can conclude that sim- 
ulation models based on RTAB tabulated S-matrix results, 
modified form factors with anomalous scattering factors and 
EPDL97 are preferable for total cross section calculations to 
models using Storm and Israel's and XCOM databases, or 
Penelope 2001 and Brennan and Cowan's parameterizations. 
It appears also justified to conclude that the cross sections 
calculated by Geant4 G4XrayRayleighModel do not represent 
a realistic physics model. 

D. Nuclear Thomson scattering 

The generation of the elastically scattered photon in gen- 
eral purpose Monte Carlo codes only accounts for Rayleigh 
scattering; RTAB tabulations enable accounting also for nu- 
clear Thomson scattering in the angular distribution of the 
scattered photon. The contribution due to the nuclear Thomson 
amplitude is more relevant at higher energies and large angles; 
as it can be observed in Fig. |2j the inclusion of nuclear 



20 



Z=13 




0.1 1 10 100 1000 



Energy (keV) 
Z=50 




0.1 1 10 100 1000 



Z=29 

600 r-i ' ' 




0.1 1 10 100 1000 



Energy (keV) 



Z=82 




0.1 1 10 100 1000 



Energy (keV) Energy (keV) 

Fig. 10. Total cross section as a function of energy for aluminium, copper, tin and lead: calculations based on S-matrix (SM, black filled circles), EPDL 
(red dashed line), modified form factors with anomalous scattering factors MFASF (blue solid line), XCOM (green dotted line), Storm and Israel (brown 
dotted-dashed line), Penelope 2001 (turquoise double-dashed line), Brennan and Cowan (yellow dashed line) and G4XrayRayleighModel (G4std, magenta 
long-dashed line). 
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Fig. 11. Total cross section as a function of energy for barium, tungsten, lead and bismuth: calculations based on S-matrix (SM, black solid line), EPDL (red 
dashed line), modified form factors with anomalous scattering factors MFASF (blue double-dashed line), XCOM (green dotted line), Storm and Israel (brown 
dotted-dashed line), Penelope 2001 (turquoise double-dashed line), Brennan and Cowan (yellow dashed line) and G4XrayRayleighModel (G4std, magenta 
long-dashed line), and experimental data (black circles) from [159] . In some plots the experimental errors are smaller than the symbol size. 
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Fig. 12. Total cross section as a function of energy for silver, cadmium, tin, iodine, cerium, gadolinium, thulium, tantalum and mercury: calculations based 
on S-matrix (SM, black solid line), EPDL (red dashed line), modified form factors with anomalous scattering factors MFASF (blue double-dashed line), 
XCOM (green dotted line) and Storm and Israel (brown dotted-dashed line), Brennan and Cowan (yellow dashed line) and G4XrayRayleighScatteringModel 
(G4std, magenta long-dashed line) and experimental measurements from [160] (black dots). The Penelope 2001 model, not shown for better clarity of the 
plots, exhibits similar behaviour to XCOM cross sections. In some plots the experimental errors are smaller than the symbol size. 
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Thomson scattering in the simulation has significant effects 
on the correct estimate of photon backscattering. 

The effect of accounting for the nuclear Thomson amplitude 
on the accuracy of the simulation is summarized in Table 



XIII The results concern the comparison of differential cross 
sections with experiment, and consist of the efficiency for 
two variants of the S-matrix model, respectively accounting 
for Rayleigh scattering and Nuclear Thomson scattering, or 
neglecting the latter; they are reported in three test con- 
figurations, concerning the whole data sample and the data 
respectively below and above 90°. In all the configurations 
the inclusion of nuclear Thomson scattering contributes to 
improving the accuracy of the simulation. 

TABLE XIII 

Efficiencies for differential cross sections based on S-matrix 
calculations accounting for rayleigh scattering only or 
including also nuclear thomson scattering 





SM 


SM 




Rayleigh 


Rayleigh + Nuclear Thomson 


all 


0.72 ± 0.05 


0.77 ± 0.05 


< 90° 


0.76 ± 0.05 


0.82 ± 0.05 


> 90° 


0.35 ± 0.12 


0.59 ± 0.12 



E. Simulation at the electronvolt scale 

Microdosimetry and nanodosimetry simulation requires the 
capability of modeling particle interactions with matter down 
to the scale of a few electronvolts. EPDL97 includes tabula- 
tions of Rayleigh total cross sections down to 1 eV, never- 
theless, its documentation warns about using the data below 
100 eV for photon transport calculations, since "the uncer- 
tainty in the data rapidly increases with decreasing energy" 
[19]. 

A qualitative comparison of EPDL97 total Rayleigh scat- 
tering cross sections and experimental data [177]-[190] for 
rare gases, hydrogen and atomic nitrogen at energies between 
1.8 and 11.2 eV is illustrated in Fig. [13] Although one 



can observe some resemblance between EPDL97 calculations 
and experimental cross section measurements in some of the 
plots and large discrepancies in other ones, a quantitative 
appraisal of the accuracy of EPDL97 is not realistic due 
to the scarcity of experimental data, their inconsistencies 
and possible systematic effects introduced by normalization 
procedures applied to some of the experimental data. 

Further experimental measurements would be needed for a 
thorough assessment of simulation models meant for micro- 
dosimetry simulation. 

F. Interpolation algorithms 

Simulation models based on data tabulations calculate the 
physical quantities needed in the course of particle transport 
by interpolation over the tabulated values. This method has the 
advantage of avoiding the computation of complex analytical 
formulae in the course of simulation execution; nevertheless, 
interpolation methods themselves could be computationally 
expensive, especially if they involve calls to library functions 



(e.g. the calculation of logarithms) rather than elementary 
arithmetic operations only. 

EPDL documentation recommends double-logarithmic in- 
terpolation over the tabulated data; this prescription has been 
adopted for the related simulation models [45]. The same 
interpolation method has been used also for the other models 
based on form factor tabulations to ensure a uniform treatment. 

The use of S-matrix tabulations requires two-dimensional 
interpolation over energy and angle for any given target atomic 
number. Three interpolation algorithms have been evaluated, 
considering both the resulting physics accuracy and their com- 
putational performance: linear and logarithmic interpolation, 
and a simplified interpolation method described in [22], which 
utilizes scaling factors related to MFASF tabulations. Table 



XIV summarizes the efficiency at producing results consistent 
with experiment associated with various interpolation method; 
they concern the Rayleigh scattering amplitude only, but the 
conclusions hold also for other amplitudes. The logarithmic in- 
terpolation algorithm determines differential cross sections that 
are compatible with experiment at 0.01 level of significance 
in a larger number of test cases; nevertheless, the efficiencies 
resulting from the three algorithms are compatible within one 
standard deviation. 

Although less accurate than logarithmic interpolation, linear 
interpolation of S-matrix tabulations produces estimates in bet- 
ter agreement with experiment than models based on the form 
factor approximation; therefore, if computational performance 
is a concern, one can opt for linearly interpolating S-matrix 
tabulations, still preserving the superior accuracy of this model 
with respect to other physics approximations despite some 
degradation with respect to logarithmic interpolation. 

Programming techniques for performance optimization 
[191]— [193] and more refined interpolation methods than those 
discussed in this paper can be applied for efficient tabulated 
data management; their in-depth discussion is outside the 
scope of this paper. 

TABLE XIV 
Effect of interpolation algorithms 





SM 


SM 


SM 




linear 


logarithmic 


simplified 


Test cases 


71 


71 


71 


Pass 


49 


51 


48 


Fail 


22 


20 


23 


Efficiency 


0.69 


0.72 


0.68 


Error 


±0.05 


±0.05 


±0.05 



G. Computational performance 

Computational performance is an important issue in as- 
sessing whether a physics model is suitable for Monte Carlo 
particle transport, and in optimizing the selection of physics 
models for the simulation of a given experimental scenario, 
when multiple options are available. 



The validation analysis documented in section V-B has 
identified the simulation model based on RTAB tabulations of 
second order S-matrix calculations as the most accurate at re- 
producing experimental measurements of elastic scattering dif- 
ferential cross sections, i.e. at determining the actual outcome 



24 




Z=18 



Z=54 




Energy (eV) 



Energy (eV) 



Z=1 



Z=7 



Z=36 





9 10 

Energy (eV) 



Fig. 13. Total cross section as a function of energy for rare gases (helium, neon, argon, krypton and xenon), hydrogen and nitrogen in the electronvolt 
energy range: EPDL (black solid line), and experimental measurements [177]-[190]. The legend reported in the plot for neon is pertinent to all the plots; the 
experimental data are identified by the name of the first author of the corresponding publication. Please note that the plots for helium, argon and xenon are 
in logarithmic scale on the vertical axis, while the others are in linear scale. The error bars that are not visible in the plots are smaller than the symbol size. 
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of the scattering process. Nevertheless, the use of S-matrix tab- 
ulations in particle transport requires more complex operations 
than models based on the form factor approximation, since it 
implies the management of two-dimensional tabulations (over 
a grid of energies and angles) and the development of an 
ad hoc algorithm to sample the direction of the scattered 
photon, replacing the sampling algorithms that are currently 
established to deal with the form factor approximation. The 
comparative evaluation of its computational performance with 
respect to models based on the form factor approximation is 
relevant to determine whether this modeling approach, that has 
not yet been exploited in Monte Carlo codes, is practically 
usable in that context. 

From another standpoint, less accurate models, for instance 
based on parameterizations or neglecting some physical fea- 
tures, could be justified for inclusion in a simulation system, 
if they provide significantly faster computational performance 
than more accurate alternatives. 

The analysis addresses various issues related to how fea- 
tures specific to different photon elastic scattering modeling 
approaches and software algorithms affect the computational 
performance of the simulation. Evaluations of computational 
performance are necessarily specific to a given Monte Carlo 
code, computing platform and test case; the results reported 
here concern test cases using the Geant4 toolkit, and have been 
obtained on an AMD 2.4 GHz 2-core processor computer with 
Scientific Linux 5 operating system and the gcc 4.1.2 compiler. 

From a computational perspective, one can distinguish two 
methods of total cross section calculation as a function of 
energy: either by interpolation of tabulated values, or by 
means of analytical formulae. In the context of Geant4 the 
former method is used by the models based on EPDL97, 
while the latter is used by G4XrayRayleighModel and by the 
model derived from Penelope 2001. Total cross sections based 
on S-matrix calculations are computed by interpolation of 
tabulations as well; therefore, in computational terms, they 
do not differ from calculations currently based on EPDL97. 

Two test cases have been evaluated as representative of 
the computational burden associated with either method; they 
concern the average time to calculate one million total cross 
sections, for atomic number between 1 and 99, and photon 
energy between 250 eV and 50 keV. The computation time is 
0.277 ±0.019 seconds for a model based on the interpolation 
of EPDL97, and 0.332 ± 0.007 seconds for the analytical 
algorithm implemented in G4XrayRayleighModel\ it appears 
independent from the photon energy. The hypothesis of equiv- 
alent computational impact of the two calculation methods 
is rejected 99% confidence level in this test case, i.e. the 
calculation of total cross sections based on the interpolation of 
EPDL97 is significantly faster that the calculation performed 
by the simplified G4XrayRayleighModel. However, this result 
should not be generalized to the computational equivalence of 
interpolation and analytical algorithms in absolute terms, since 
the relative performance of these calculation methods depends 
on the computational complexity of the analytically formulae 
and on the characteristics of the interpolation algorithms 
subject to test. 

Nevertheless, the computational impact of total cross section 



calculation methods on Monte Carlo simulation is marginal, 
since in general it affects the computational performance only 
in the initialization phase of the simulation: for instance, in 
Geant4 atomic total cross sections are calculated at initial- 
ization in the process of building mean free path tabulations 
for each material present in the geometrical model of the 
experimental set-up, which are then used in the course of 
particle transport. 

Two issues contribute to the computational performance 
of the actual photon scattering simulation: the calculation of 
the physical distribution to sample the scattering angle, and 
the sampling algorithm. The calculation of the distribution 
to sample can occur by interpolation of tabulated values or 
through an analytical formula. 

The computational performance of a set of representative 
methods for the generation of the elastically scattered photon 
has been analyzed on a test case consisting of one million pri- 
mary photons with energy between 250 eV and 1 MeV, which 
are only subject to elastic scattering as possible interaction 
with matter. The Geant4-based application developed for this 
test can be configured to use any of the existing Rayleigh scat- 
tering options in Geant4 9.5 or the newly developed software 
described in this paper, namely elastic scattering based on S- 
matrix calculations; it provides options for elemental targets 
with atomic number between 1 and 99. 



The results of this test are illustrated in Fig. 14 which report 
the average time to simulate one million events as a function 
of photon energy for two representative test cases (carbon and 
lead as targets); the plotted values are scaled to the time needed 
to simulate one million events with 10 keV primary photons 
and copper as a target. 

The computational perfomance of physics models based on 
the form factor approximation appears independent from the 
way the sampling distribution is calculated, but is strongly 
affected by the efficiency of the sampling algorithm. It is 
comparable for the Penelope 2001 and 2008 scattering imple- 
mentations, which adopt the same sampling algorithm (based 
on inverse transform adaptive sampling [48]), but calculate 
the sampling distribution respectively from an analytical for- 
mula and by interpolation. It exhibits large differences for 
the Geant4 Penelope 2008 and Livermore implementations, 
which both interpolate EPDL97 tabulations, but use different 
sampling algorithms (acceptance-rejection sampling for the 
Livermore model, a form of inverse transform sampling for 
Penelope [40]). The inefficiency of the acceptance-rejection 
algorithm shows a strong dependence on the photon energy, 
which severely affects the computational performance of the 
elastic scattering simulation above approximately a few tens 
keV. 

The simplified physics model implemented in 
G4XrayRayleighModel does not exhibit any significant 
computational advantage with respect to the interpolation of 
EPDL97 with inverse transform sampling, which is available 
in Geant4 through the reengineered Penelope 2008 model. 

The generation of photon elastic scattering based on S- 
matrix calculations is approximately a factor two slower than 
the method based on the form factor approximation exploiting 
inverse transform sampling, which, among the test cases illus- 
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trated in Figs. 14 is represented by the reengineered Penelope 
2008 model: the difference is due to performing interpolation 
over bi-dimensional tabulations rather than one-dimensional 
interpolation. The final state generation model based on S- 
matrix calculations exploits the inverse transform algorithm 
to sample the scattering angle, similarly to the Penelope 
models; thanks to its more efficient sampling algorithm, it is 
significantly faster than the model based on EPDL97 form 
factors currently implemented in Geant4, which utilizes an 
acceptance-rejection sampling algorithm. 

A factor two degradation in the computational performance 
of photon elastic scattering simulation can be considered an 
acceptable compromise for most experimental applications 
requiring accurate simulation of photon interactions, given the 
superior physical performance of this model demonstrated by 
the validation documented in section IV-BI 

It should be stressed that the degradation of performance 
discussed in this context is limited to photon elastic scattering, 
while in experimental practice other types of photon interac- 
tions, and of their secondary products, are usually activated 
in the simulation: the factor two penalty in computational 
performance represents an upper limit, which corresponds 
to applications only involving photon elastic scattering. The 
relative impact of photon elastic scattering on the overall 
computational performance cannot be quantified in general 
terms, since it depends on many physical and user parameters 
specific to a given experimental scenario: some factors affect- 
ing the computational speed of a simulation are the energy 
of the particles involved, the materials of the experimental 
set-up, the physics processes and models activated in the 
simulation, secondary production thresholds and other user 
cuts, the complexity of the geometrical configuration, and 
the characteristics of the user application software, but this 
enumeration is not intended to be exhaustive. 

A policy class for final state generation has been developed 
in the context of this study, which can exploit any form 
factor tabulation and uses an inverse transform algorithm to 
sample the photon scattering angle. This class provides optimal 
computational performance, while it substantially improves 
simulation accuracy with respect to the Penelope 2008 model, 
still in the form factor approximation, by using modified form 
factors with angle independent anomalous scattering factors. 

VI. Conclusion 

An extensive set of models for the simulation of photon 
elastic scattering has been quantitatively evaluated regarding 
their accuracy at reproducing experimental measurements and 
their computational performance. 

The analysis has identified the simulation of photon elas- 
tic scattering based on second order S -matrix calculations 
tabulated in RTAB as the state-of-the-art. This model, that 
accounts for Rayleigh scattering and nuclear Thomson scat- 
tering, represents approximately a factor two improvement in 
compatibility with experiment with respect to the Rayleigh 
scattering models currently available in Geant4 and in other 
general purpose Monte Carlo codes. The inclusion of nuclear 
Thomson scattering is relevant especially at higher energies 



and for backscattering. Complementary evaluation of the com- 
putational performance has demonstrated that photon elastic 
scattering simulation based on S -matrix calculations is practi- 
cally feasible in a general purpose Monte Carlo system, despite 
its greater computational complexity than models based on the 
form factor approximation currently in use. 

The results of the analysis hint that the accuracy of simula- 
tion based on S -matrix calculations could be further improved 
by optimizing the tabulation grid of RTAB for more precise in- 
terpolation in areas that are sensitive to the underlying atomic 
structure. Nevertheless, the production of an extended data 
library of second order S -matrix calculations would require 
significant computational investment and the collaboration of 
expert theoreticians. 

If one prefers to base the simulation of Rayleigh scat- 
tering on the form factor approach to avoid more complex 
computations required by the S -matrix model, modified form 
factors with anomalous scattering factors (MFASF) appear 
the preferable choice among the various form factor options 
examined in this paper, although they result in degraded 
simulation accuracy with respect to the state-of-the-art model 
based on S -matrix calculations. The choice of the method 
to sample the direction of the scattered photon is critical 
for computational performance: inverse transform sampling is 
preferable to acceptance-rejection. 

Relativistic form factors result in worse accuracy than other 
form factor options (non-relativistic and modified); therefore 
their use in simulation models is not encouraged. 

Anomalous scattering factors calculated in the RTAB theo- 
retical environment contribute to improve the accuracy of cross 
sections based on modified and relativistic form factors, while 
the anomalous scattering factors included in EPDL97 do not 
improve the results based on this data library. 

Total and differential cross sections calculated by interpola- 
tion of pristine EPDL97 tabulations, and by the reengineered 
code and reprocessed EPDL97 tabulations originating from 
Penelope 2008 produce equivalent physics results. The main- 
tenance in Geant4 of two implementations that conform to the 
same software design scheme and produce equivalent physics 
outcome appears redundant. 

The total cross sections calculated by the Penelope 2001 
model neglect the underlying atomic structure in the low 
energy range, and the differential cross sections exhibit signif- 
icantly worse accuracy than other models based on the form 
factor approximation: the maintenance of the Penelope 2001 
model in Geant4 does not appear physically justified. 

Total cross sections calculated by XCOM, Storm and Israel, 
and Brennan and Cowan look qualitatively similar; they exhibit 
relatively limited differences with respect to cross sections 
based on EPDL97 and on S -matrix calculations above approxi- 
mately 10 keV, but they appear largely deficient at reproducing 
effects close to absorption edges. 

The total cross sections calculated by Geant4 
G4XrayRayleighModel are insensitive to the atomic structure 
at low energies, and largely inconsistent with experimental 
data in the energy range of a few hundred keV. The scattering 
angle distribution produced by this model corresponds 
to Thomson scattering by a point-like charge, which is an 
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Fig. 14. Time spent to simulate one million photon elastic scatterings with carbon (left) and lead (right), as a function of photon energy with current Geant4 
Rayleigh scattering implementations and with the new model based on S-matrix calculations: reengineered Penelope 2001 and 2008 models (respectively green 
diamond and blue cross), EPDL-based "Livermore" model (red circles), G4XrayRayleighModel (magenta asterisk) and S-matrix model (SM, black squares). 
The timing is expressed in relative units with respect to the time spent to simulate one million scatterings with copper at 10 keV using the Geant4 model 
reengineered from Penelope 2008. The error bars are not visible in the plot, as they are smaller than the symbol size. 



unrealistic description of photon elastic scattering. This model 
does not exhibit any computational advantage with respect to 
more accurate models currently available in Geant4. Further 
maintenance of this model in Geant4 should be seriously 
considered. 

The policy-based software design adopted in this study for 
the simulation of photon elastic scattering has played a key 
role in ensuring versatility of modeling, at the same time 
greatly facilitating the verification and validation process due 
to minimization of the dependencies of basic physics code 
on other parts of Geant4. This finding is relevant in view 
of future design evolutions of the Geant4 toolkit to improve 
its reliability and to reduce the efforts invested for quality 
assurance and maintenance. 

The software for the simulation of photon elastic scattering 
developed for this study will be proposed for release in a 
forthcoming Geant4 version, following the publication of this 
paper, to improve Geant4 capabilities in this physics domain. 

Acknowledgment 

The authors are grateful to Lynn Kissel for valuable infor- 
mation and advice, and to Andrew Buckley for proofreading 
the manuscript and useful comments. 

The CERN Library, in particular Tullio Basaglia, has pro- 
vided helpful assistance and essential reference material for 
this study. 

References 

[1] R. Cesareo, A. L. Hanson, G. E. Gigante, L. J. Pedraza, and S. Q. 
G. Mahtaboally, "Interaction of keV photons with matter and new 
applications", Phys. Rep., vol. 213, no. 3, pp. 117-178, 1992. 



[2] J. H. Hubbell, "Review and history of photon cross section calcula- 
tions", Phys. Med. Biol, vol. 51, pp. R245-R262, 2006. 

[3] P. P. Kane, L. Kissel, R. H. Pratt, and S. C. Roy, "Elastic scattering of 
7-rays and X-rays by atoms", Phys. Rep., vol. 140, no. 2, pp. 75-159, 
1986. 

[4] S. C. Roy, L. Kissel, and R. H. Pratt, "Elastic scattering of photons", 

Rad. Phys. Chem. , vol. 56, pp. 3-26, 1999. 
[5] D. A. Bradley, O. D. Goncalves, and P. P. Kane, "Measurements of 

photon-atom elastic scattering cross-sections in the photon energy range 

1 keV to 4 MeV", Rad. Phys. Chem. , vol. 56, pp. 125-150, 1999. 
[6] K. Amako et al., "Comparison of Geant4 electromagnetic physics 

models against the NIST reference data", IEEE Trans. Nucl. Set, vol. 

52, no. 4, pp. 910-918, 2005. 
[7] T. G. Trucano, L. P. Swiler, T. Igusa, W. L. Oberkampf, and M. Pilch, 

"Calibration, validation, and sensitivity analysis: Whats what", Reliab. 

Eng. Syst. Safety, vol. 91, no. 10-11, pp. 1331-1357, 2006. 
[8] G. A. P. Cirrone, G. Cuttone, F. Di Rosa, L. Pandola, F. Romano, and 

Q. Zhang, "Validation of the Geant4 electromagnetic photon cross- 
sections for elements and compounds", Nucl. Instr. Meth. A, vol. 618, 

pp. 315-322, 2010. 
[9] S. Agostinelli et al., "Geant4 - a simulation toolkit" Nucl. Instrum. 

Meth. A, vol. 506, no. 3, pp. 250-303, 2003. 
[10] J. Allison et al., "Geant4 Developments and Applications" IEEE Trans. 

Nucl. Sci., vol. 53, no. 1, pp. 270-278, 2006. 
[11] R. H. Pratt, L. Kissel, and P. M. Bergstrom, "New relativistic S- 

matrix results for scattering beyond the usual anomalous factors/beyond 

impulse approximation", in Resonant Anomalous X-ray Scattering - 

Theory and Applications, Elsevier Science, Amsterdam, pp. 9-33, 1994. 
[12] N. A. Dyson, "X-rays in Atomic and Nuclear Physics", Cambridge 

University Press, 2nd ed., pp. 215-219, 1990. 
[13] J. H. Hubbell, W. J. Viegele, E. A. Briggs, R. T. Brown, D. T. Cromer, 

and R. J. Howerton, . "Atomic form factors, incoherent scattering 

functions, and photon scattering cross sections", /. Phys. Chem. Ref. 

Data, vol. 4, pp. 471-538, 1975. 
[14] J. H. Hubbell and I. Overb0, "Relativistic atomic form factors and 

photon coherent scattering cross sections", /. Phys. Chem. Ref. Data, 

vol. 8, pp. 69-105, 1979. 
[15] D. Schaupp, M. Schumacher, F. Smend, P. Rullhusen, and J. H. 

Hubbell, "Small angle Rayleigh scattering of photons at high energies: 



28 



Tabulation of relativistic HFS modified atomic form factors", /. Phys. 
Chem. Ref. Data, vol. 12, pp. 467-508, 1983. 

G. E. Brown, R. E. Peierls, J. B. Woodward, "The coherent scattering 
of 7-rays by K-electrons in heavy atoms", Proc. R. Soc. London A, 
vol. 227, no. 1168, pp. 51-58, 1955. 

L. Kissel and R. H. Pratt, "Rayleigh scattering - elastic photon 
scattering by bound electrons", in Atomic Inner-shell Physics, Plenum 
Press, New York, pp. 465-498, 1985. 

L. Kissel, B. Zhou, S. C. Roy, S. K. Sen Gupta, and R. H. Pratt, "The 
validity of form factor, modified form factor and anomalous scattering 
factor approximations in elastic scattering calculations" Acta Cryst. A, 
vol. 51, pp. 271-288, 1995. 

D. Cullen et al, "EPDL97, the Evaluated Photon Data Library", 
Lawrence Livermore National Laboratory Report UCRL-50400, Vol. 

6, Rev. 5, 1997. 

M. B. Chadwick et al., "ENDF/B-VII.l Nuclear Data for Science and 
Technology: Cross Sections, Covariances, Fission Product Yields and 
Decay Data", Nucl. Data Sheets, vol. 112, no. 12, pp. 2887-3152, 2011. 

D. E. Cullen, "Program SCATMAN: A Code Designed to Calculate 
Photon Coherent Scattering Anomalous Scattering Factors and Cross 
Sections", Lawrence Livermore National Laboratory UCRL-ID- 103422 
Report, 1989. 

L. Kissel, "RTAB: the Rayleigh scattering database", Radiat. Phys. 
Chem., vol. 59, pp. 185-200, 2000. 

L. R. M. Morin, "Molecular Form Factors and Photon Coherent 
Scattering Cross Sections of Water", /. Phys. Chem. Ref. Data, vol. 
11, pp. 1091-1098, 1982. 

H. Hirayama, Y. Namito, A. F. Bielajew, S. J. Wilderman, and W. 
R. Nelson, "The EGS5 Code System", SLAC-R-730 Report, Stanford, 
CA, 2006. 

W. R. Nelson, H. Hirayama, and D. W. O. Rogers, "The EGS4 Code 
System", SLAC-265 Report, Stanford, CA, 1985. 

E. Storm and H. I. Israel, "Photon cross sections from 1 keV to 100 
MeV for elements Z=l to Z=100", Atom. Data Nucl. Data Tables, vol. 

7, pp. 565-681, 1970. 

H. P. Hanson, F. Herman, J. D. Lea, and S. Skillman, "HFS Atomic 
Scattering Factors", Acta Cryst., vol. 17, pp. 1040-1044, 1964. 

I. Kawrakow ,E. Mainegra-Hing, D.W.O. Rogers, F. Tessier and B.R.B. 
Walters, "The EGSnrc Code System: Monte Carlo Simulation of 
Electron and Photon Transport NRCC PIRS-701, 5th printing, 2010. 
M. J. Berger et al., "XCOM: Photon Cross section Database (version 
1.5)", National Institute of Standards and Technology, Gaithersburg, 
MD, 2010. [Online] Available: http://physics.nist.gov/xcom 

S. M. Seltzer, "Electron-Photon Monte Carlo Calculations: The 
ETRAN Code", Appl. Radiat. Isot., vol. 42, no. 10, pp. 917-941, 1991. 
G. Battistoni et al., "The FLUKA code: description and benchmarking", 
AIP Conf. Proc, vol. 896, pp. 31-49, 2007. 

A. Ferrari et al., "Fluka: a multi-particle transport code", Report CERN- 
2005-010, INFN/TC-05/11, SLAC-R-773, Geneva, Oct. 2005. 

B. C. Franke, R. P. Kensek and T. W. Laub, "ITS5 theory manual", rev. 
1.2, Sandia Natl. Lab. Report SAND2004-4782, Albuquerque, 2004. 
X-5 Monte Carlo Team, "MCNP - A General Monte Carlo N-Particle 
Transport Code, Version 5", Los Alamos National Laboratory Report 
LA-UR-03-1987, 2003, revised 2008. 

J. S. Hendricks et al., "MCNPX, Version 2.6.0", Los Alamos National 
Laboratory Report LA-UR-06-7991, 2008. 

D. E. Cullen et al., "Tables and Graphs of Photon Interaction Cross 
Sections from 10 eV to 100 GeV Derived from the LLNL Evaluated 
Photon Data Library (EPDL)", Lawrence Livermore National Labora- 
tory Report UCRL-50400, Vol. 6, Rev. 4, 1989. 
D. Garber Ed. , "ENDF/B Summary Documentation", BNL-17541, 
1975. 

D. B. Pelowitz et al., "MCNPX 2.7.E Extensions", Los Alamos 
National Laboratory Report LA-UR-1 1-0 1502, 2011. 
J. Baro, J. Sempau, J. M. Fernandez- Varea, and F. Salvat, "PENELOPE, 
an algorithm for Monte Carlo simulation of the penetration and energy 
loss of electrons and positrons in matter", Nucl. Instrum. Meth. B, vol. 
100, no. 1, pp. 31-46, 1995. 

F. Salvat, J.M. Fernandez- Varea, and J. Sempau, "Penelope - A code 
system for Monte Carlo simulation of electron and photon transport", 
Proc. Workshop NEA 6416, 2008. 

F. Salvat, J.M. Fernandez-Varea, and J. Sempau, "Penelope-2011 - 
A code system for Monte Carlo simulation of electron and photon 
transport", Proc. Workhop NEA/NSC/DOC(2011)5, 2011. 
J. Baro, M. Roteta, J. M. Fernandez-Varea, and F. Salvat, "Analytical 
cross sections for Monte Carlo simulation of photon transport", Radiat. 
Phys. Chem., vol. 44, no. 5, pp. 531-552, 1994. 



[43 

[44; 

[45 
[46 

[47 
[48 

[49 
[50 



[51 

[52 



[53 
[54 

[55 

[56 
[57 

[58 

[59 
[60 

[61 

[62 
[63 
[64 
[65 
[66 
[67 
[68 

[69 



"GEANT Detector Description and Simulation Tool", CERN Program 
Library Long Writeup W5013, 1995. 

R. L. Ford and W. R. Nelson, "The EGS code system: computer 
programs for the Monte Carlo simulation of electromagnetic cascade 
showers (version 3)", SLAC-210 Report, Stanford, CA, 1978. 
J. Apostolakis, S. Giani, M. Maire, P. Nieminen, M.G. Pia, L. Urban, 
"Geant4 low energy electromagnetic models for electrons and photons" 
INFN/AE-99/18, Frascati, 1999. 

S. Chauvie, G. Depaola, V. Ivanchenko, F. Longo, P. Nieminen and 
M. G. Pia, "Geant4 Low Energy Electromagnetic Physics", in Proc. 
Computing in High Energy and Nuclear Physics, Beijing, China, pp. 
337-340, 2001. 

S. Chauvie et al., "Geant4 Low Energy Electromagnetic Physics", in 
2004 IEEE Nucl. Sci. Symp. Conf. Rec, pp. 1881-1885, 2004. 

F. Salvat, J.M. Fernandez-Varea, E. Acosta, and J. Sempau, "Penelope 
- A code system for Monte Carlo simulation of electron and photon 
transport", Proc. Workshop NEA, 2001. 

S. Brennan and P. L. Cowan, "A suite of programs for calculating 
xray absorption, reflection, and diffraction performance for a variety 
of materials at arbitrary wavelengths", Rev. Sci. Instrum., vol., 63, pp. 
850-853, 1992. 

W. H. McMaster, N. Kerr Del Grande, J. H. Mallet and J. H. 
Hubbell, "Compilation of X-ray cross sections", Section II Revision 1,, 
Lawrence Livermore National Laboratory Report UCRL-50174, 1969. 
A. Alexandrescu, "Modern C++ Design", Ed.: Addison- Wesley, 2001. 
S. Chauvie et al., "Geant4 physics processes for microdosimetry 
simulation: design foundation and implementation of the first set of 
models", IEEE Trans. Nucl. Sci., vol. 54, no. 6, pp. 2619-2628, 2007. 
R. Y. Rubinstein and D. P. Kroese, "Simulation and the Monte Carlo 
Method", 2 nd ed., Wiley, Hoboken, NJ, 2008. 

M. Augelli et al., "Research in Geant4 electromagnetic physics design, 
and its effects on computational performance and quality assurance", 
in 2009 IEEE Nucl. Sci. Symp. Conf. Rec, pp. 177-180, 2009. 
M. G. Pia et al., "Design and performance evaluations of generic 
programming techniques in a R&D prototype of Geant4 physics", J. 
Phys.: Conf. Ser, vol. 219, pp. 042019, 2010. 

M. Mitchell, Engauge Digitizer [Online]. Available: 
http://digitizer.sourceforge.net 

A. Wald and J. Wolfowitz, "An exact test for randomness in the non- 
parametric case, based on serial correlation," Ann. Math. Stat., vol. 14, 
pp. 378-388, 1943. 

L. W. Alvarez, F. S. Crawford, and M. L. Stevenson, "Elastic scattering 
of 1.6 MeV gamma rays from H, Li, C, and Al nuclei", Phys. Rev., 
vol. 112, no. 4, pp. 1267-1273, 1958. 

S. Anand, M. Singh, and B. S. Sood, "Large angle Rayleigh scattering 
of 662 keV gamma rays", Curr. Sci., vol. 32, no. 9, pp. 401-402, 1963. 
S. Anand and B. S. Sood, "Measurement of Z dependence of elastic 
scattering cross section of gamma rays", Nucl. Phys., vol. 73, no. 2, 
pp. 368-378, 1965. 

J. Banaigs, P. Eberhard, L. Goldzahl, and E. Hara, "Diffusion elastique 
des rayons 1.12-1.17-1.33 et 2.62 MeV", /. Phys. Radium, vol. 19, no. 
1, pp. 70-72, 1958. 

C. Baraldi, E. Casnati, A. Tartari, M. Andreis, and B. Singh, "Mea- 
surements of anomalous elastic scattering of 59.54 keV photons", Phys. 
Rev. A, vol. 54, no. 6, pp. 4947-4953, 1996. 

S. de Barros, O. Goncalves, M. Gaspar, and J. R. Moreira, "Elastic 
scattering of 412, 468, and 662 keV 7 rays", Phys. Rev. C, vol. 22, 
no. 1, pp. 332-336, 1980. 

S. de Barros, J. Eichler, M. Gaspar, and O. Goncalves, "Rayleigh 
scattering of 468 keV photons by different atoms", Phys. Rev. C, vol. 
24, no. 4, pp. 1765-1768, 1981. 

S. de Barros, J. Eichler, O. Gonccalves, and M. Gaspar, "Experiments 
on Rayleigh scattering of 145 keV and 317 keV photons", Z. Natur- 
forsch. A, vol. 36, pp. 595-599, 1981. 

G. Basavaraju, P. Kane, L. Kissel, and R. Pratt, "Elastic scattering of 81 
keV 7 rays by aluminium, nickel, tantalum, gold and lead", Pramana, 
vol. 44, pp. 545-553, 1995. 

G. Basavaraju, P. Kane, and K. Varier. "Elastic scattering of 1.17 and 
1.33 MeV gamma rays by molybdenum, tantalum and lead", Pramana, 
vol. 12, pp. 665-678, 1979. 

G. Basavaraju and P. P. Kane, "Elastic scattering of 1.12 MeV and 
1.33 MeV 7-rays through 90° and 124.5° by different elements and 
the real parts of the Delbr'uck scattering amplitudes", Nucl. Phys. A, 
vol. 149, no. 1, pp. 49-62, 1970. 

G. Basavaraju, P. P. Kane, L. D. Kissel, and R. H. Pratt, "Elastic 



29 



scattering of 81 keV 7 rays", Phys. Rev. A, vol. 49, no. 5, pp. 3664- 
3672, 1994. 

[70] G. Basavaraju, P. P. Kane, S. M. Lad, L. Kissel, and R. H. Pratt, "Elastic 
scattering of 88.03 keV 7 rays", Phys. Rev. A, vol. 51, no. 3, pp. 2608- 
2610, 1995. 

[71] A. M. Bernstein and A. K. Mann, "Scattering of gamma rays by a 
static electric field", Phys. Rev., vol. 110, no. 4, pp. 805-814, 1958. 

[72] D. A. Bradley, C. S. Chong, A. A. Tajuddin, A. Shukri, and A. M. 
Ghose, "Small-angle coherent gamma-ray scattering at 662 keV for Pb 
and Sn", Nucl. Instr. Meth. A, vol. 280, no. 2-3, pp. 207-211, 1989. 

[73] D. A. Bradley, C. S. Chong, A. A. Tajuddin, A. Shukri, and A. M. 
Ghose, "Small-angle coherent 7-ray scattering for moderate- to high- 
atomic-number elements", Phys. Rev. A, vol. 41, no. 11, pp. 5974-5979, 
1990. 

[74] D. A. Bradley and A. M. Ghose, "Total-atom differential coherent- 
scattering cross-section measurements on Sn and Pb using moderate- 
energy 7 rays", Phys. Rev. A, vol. 33, no. 1, pp. 191-204, 1986. 

[75] C. Bui and M. Milazzo, "Measurements of anomalous dispersion in 
Rayleigh scattering of characteristic X-ray fluorescence", Nuovo Cim. 
D, vol. 11, pp. 655-686, 1989. 

[76] E. Casnati, C. Baraldi, and A. Tartari, "Measurements of the total 
atomic differential cross section of elastic scattering of 59.54 keV 
photons", Phys. Rev. A, vol. 42, no. 5, pp. 2627-2633, 1990. 

[77] B. K. Chatterjee and S. C. Roy. "Tables of elastic scattering cross 
sections of photons in the energy range 50-1500 keV for all elements 
in the range 13 ^ Z ^ 104", /. Phys. Chem. Ref. Data, vol. 27, no. 6, 
pp. 1011-1215, 1998. 

[78] W. Chitwattanagorn, S. Kwaengsobha, J. Jarasrangsichol, K. Rach- 
dawanpong, and P. Oungmanee, "Elastic scattering of 152 Eu 7 rays 
by various targets", Nucl. Instr. Meth. A, vol. 255, no. 1-2, pp. 75-77, 
1987. 

[79] W. Chitwattanagorn, R. B. Taylor, P. Teansomprasong, and I. B. 
Whittingham, "Elastic scattering of 152 Eu 7rays by Pb", /. Phys. G, 
vol. 6, no. 9, pp. 1147-1165, 1980. 

[80] C. S. Chong, A. A. Tajuddin, A. Shukri, and D. A. Bradley, "Photon- 
atom scattering of 26.4 and 59.5 keV gamma rays in Mo, Nb and Zr", 
Radiat. Phys. Chem., vol. 47, no. 5, pp. 689-690, 1996. 

[81] N. Cindro and K. Ilakovac, "Elastic scattering of gamma-rays", Nucl. 
Phys., vol. 5, pp. 647-652, 1958. 

[82] W. R. Dixon and R. S. Storey, "The elastic scattering of 60 Co 7 rays 
from lead", Can. J. Phys., vol. 46, no. 10, pp. 1153-1161, 1968. 

[83] J. Eichler and S. de Barros, "Coherent scattering of 59.54 keV 7 rays 
by Al, Cu, Zn, Cd, and Pb", Phys. Rev. A, vol. 32, no. 2, pp. 789-792, 
1985. 

[84] J. Eichler, S. de Barros, O. Goncalves, and M. Gaspar, "Comparison 
of Compton and Rayleigh scattering at 145 keV", Phys. Rev. A, vol. 
28, no. 6, pp. 3656-3658, 1983. 

[85] I. S. Elyaseery, A. Shukri, C. S. Chong, A. A. Tajuddin, and D. A. 
Bradley, "Measurement of forward scattering of 59.5 keV gamma rays 
by low to high Z elements", Radiat. Phys. Chem., vol. 51, no. 4-6, pp. 
365-367, 1998. 

[86] I. S. Elyaseery, A. Shukri, C. S. Chong, A. A. Tajuddin, and D. A. 

Bradley, "Photon-atom scattering of 13.95-, 17.75-, 26.36-, and 59.54- 

keV photons by Cu, Zn, Zr, Nb, Mo, Ag, Cd, In, Sn, Ta, and W", 

Phys. Rev. A, vol. 57, no. 5, pp. 3469-3477, 1998. 
[87] I. S. Elyaseery, A. Shukri, C. S. Chong, A. A. Tajuddin, and D. 

A. Bradley, "Total differential scattering cross-sections of 59.54 keV 

photons for elements in the range 29 ^ Z ^ 74", Radiat. Phys. Chem. , 

vol. 59, no. 2, pp. 211-214, 2000. 
[88] S. Erzeneoglu, L. Demir, and Y. Sahin, "Experimental studies on the 

coherent scattering of 59.5 keV 7-rays by Fe, Zn and Nb", Phys. Scr., 

vol. 56, no. 1, pp. 89, 1997. 
[89] S. Erzeneoglu, Y. Kurucu, R. Durak, and Y. Sahin, "Coherent scattering 

of 59.5 keV photons by Au and Pb", Phys. Scr., vol. 54, no. 2, pp. 153, 

1996. 

[90] M. L. Garg, R. R. Garg, F. Hennrich, and D. Heimermann, "Elastic 
and inelastic scattering of photons in the X-ray energy region", Nucl. 
Instr. Meth. B, vol. 73, no. 2, pp. 109-114, 1993. 

[91] M. Gaspar, O. Goncalves, S. de Barros, and J. Eichler, "Some exper- 
iments on 7-scattering between 122 and 145 kev at small angles", Z. 
Phys. D, vol. 1, pp. 287-290, 1986. 

[92] L. Goldzahl and P. Eberhard, "Contribution du champ electrique des 
noyaux a la diffusion elastique des 7 du 60 Co", /. Phys. Radium, vol. 
18, pp. 33-43, 1957. 

[93] O. Goncalves, S. de Barros, M. Gaspar, A. M. Goncalves, and J. Eich- 
ler, "Elastic scattering of 662 kev 7-rays on atoms", Z. Phys. D, vol. 
1, pp. 167-170, 1986. 



[94] O. D. Goncalves and S. D. Magalhaes, "Accurate Rayleigh differential 

cross-sections for 60 keV photons. The Ag case", Radiat. Phys. Chem., 

vol. 59, no. 2, pp. 201-209, 2000. 
[95] O. D. Goncalves, H. Schechter, M. I. Lopes, and V. Chepel, "Rayleigh 

to Compton differential cross-section ratios in liquid xenon", X-Ray 

Spectrom., vol. 28, no. 5, pp. 384-387, 1999. 
[96] M. Gowda, S. J. Anasuya, and K. S. Puttaswamy, "Rayleigh scattering 

of 145- and 279-keV 7 rays in Al, Cu, Sn, and Pb at forward angles", 

Phys. Rev. A, vol. 33, no. 4, pp. 2301-2304, 1986. 
[97] M. L. Garg, R. R. Garg, F. Hennrich, and D. Heimermann, "Elastic 

and inelastic scattering of photons in the X-ray energy region", Nucl. 

Instr. Meth. B, vol. 73, no. 2, pp. 109114, 1993. 
[98] L. P. Guy, R. B. Taylor, and I. B. Whittingham, "Small angle elastic 

scattering of 122 and 245 keV 7-rays", Nucl. Instr. Meth. B, vol. 71, 

no. 4, pp. 361-365, 1992. 
[99] E. Hara, J. Banaigs, P. Eberhard, and J. Goldzahl, L. abd Mey, 

"Diffusion elastique par le plomb et l'etain des rayons 7 de 1.33 MeV 

et de 1.17 MeV du 60 Co, et de 1.12 MeV du 65 Zn", /. Phys. Radium, 

vol. 19, pp. 668-673, 1958. 
[100] G. Hardie, J. S. DeVries, and C. Chiang, "Elastic Scattering of 1.33 

MeV Photons from Lead and Uranium", Phys. Rev. C, vol. 3, no. 3, 

pp. 1287-1293, 1971. 
[101] G. Hardie, W. J. Merrow, and D. R. Schwandt, "Large-Angle Elastic 

Scattering of 1.33-MeV Photons from Lead and Uranium", Phys. Rev. 

C, vol. 1, no. 2, pp. 714-720, 1970. 
[102] U. Hauser and B. Mussgnug, "Elastic scattering of 0.5 MeV gamma 

rays by heavy elements", Z. Phys. A, vol. 195, pp. 252-272, 1966. 
[103] O. Icelli and S. Erzeneoglu, "Coherent scattering of 59.5 keV 7-rays 

by 79 Au through angles from 45° to 125°", Spectrochim. Acta B, vol. 

56, no. 3, pp. 331-335, 2001. 
[104] S. Kahane, R. Moreh, and O. Shahal, "Rayleigh scattering of neutron 

capture 7 rays from lead", Phys. Rev. A, vol. 46, no. 5, pp. 2489-2494, 

1992. 

[105] S. Kahane, R. Moreh, and O. Shahal, "Rayleigh scattering of E = 465- 

2842 keV neutron capture 7-rays from Ta, In and Cu", Radiat. Phys. 

Chem., vol. 48, no. 4, pp. 419-426, 1996. 
[106] S. Kahane, O. Shahal, and R. Moreh, "Rayleigh scattering of neutron 

capture 7 rays from U", Phys. Rev. A, vol. 40, no. 12, pp. 6950-6957, 

1989. 

[107] P. P. Kane, G. Basavaraju, S. M. Lad, K. M. Varier, L. Kissel, and 
R. H. Pratt, "Inelastic and anomalous elastic scattering of 88.03-keV 
7 rays", Phys. Rev. A, vol. 36, no. 12, pp. 5626-5631, 1987. 

[108] P. P. Kane, G. Basavaraju, J. Mahajani, and A. K. Priyadarsini, "Elastic 
and Compton scattering of 1.17 and 1.33 MeV gamma rays through 
small angles", Nucl. Instr. Meth., vol. 155, no. 3, pp. 467-474, 1978. 

[109] P. P. Kane, J. Mahajani, G. Basavaraju, and A. K. Priyadarsini, 
"Scattering of 1.1732- and 1.3325-MeV gamma rays through small 
angles by carbon, aluminum, copper, tin, and lead", Phys. Rev. A, vol. 
28, no. 3, pp. 1509-1516, 1983. 

[110] B. Kasten, D. Schaupp, P. Rullhusen, F. Smend, M. Schumacher, and 
L. Kissel, "Coulomb correction effect in Delbr'uck scattering and 
atomic Rayleigh scattering of 1-4 MeV photons", Phys. Rev. C, vol. 
33, no. 5, pp. 1606-1615, 1986. 

[Ill] A. Kumar, J. Shahi, D. Mehta, and N. Singh, "Differential cross-section 
measurements for elastic scattering of 22.1 keV photons by elements 
with 6 ^ Z ^ 81", Nucl. Instr. Meth. B, vol. 194, no. 2, pp. 105-111, 
2002. 

[112] A. Kumar, J. S. Shahi, M. L. Garg, S. Puri, D. Mehta, and N. Singh, 
"Large-angle elastic scattering of 88.03 keV photons by elements with 
30< Z < 92", Nucl. Instr. Meth. B, vol. 183, no. 3-4, pp. 178-188, 
2001. 

[113] S. Kumar, V. Sharma, J. S. Shahi, D. Mehta, and N. Singh, "Elastic 
scattering of 59.54-keV 7-rays in elements with 22 < Z < 92 at 
momentum transfer 0.4 < x < 4.7 A -1 ", Eur. Phys. J. D, vol. 55, 
pp. 23-33, 2009. 

[114] J. P. Lestone, R. B. Taylor, P. Teansomprasong, and I. B. Whittingham, 
"Elastic scattering of 152 Eu and 154 Eu 7 rays by atoms", Phys. Rev. 
A, vol. 37, no. 9, pp. 3218-3231, 1988. 

[115] A. C. Mandal, D. Mitra, M. Sarkar, and D. Bhattacharya, "Differential 
elastic scattering cross sections of 22.1 keV x ray by elements in the 
range 22 < Z < 82", Phys. Rev. A, vol. 66, no. 4, pp. 042705, 2002. 

[116] A. K. Mann. "Elastic scattering of gamma rays", Phys. Rev., vol. 101, 
no. 1, pp. 4-8, 1956. 

[117] S. Messelt and A. Storruste, "Elastic scattering of gamma rays", Proc. 
Phys. Soc. A, vol. 69, no. 5, pp. 381-387, 1956. 



30 



[118] W. Muckenheim and M. Schumacher, "Delbr'uck and Rayleigh scat- 
tering by uranium investigated at photon energies between 0.1 and 1.5 
MeV", /. Phys. G, vol. 6, no. 10, pp. 1237-1251, 1980. 

[119] V. A. N. Murty, V. Lakshminarayana, and S. Jnanananda, "Elastic 
scattering of 1.12 MeV gamma rays", Nucl. Phys., vol. 62, no. 2, pp. 
296-304, 1965. 

[120] V. A. N. Murty, V. Lakshminarayana, and S. Jnanananda, "The elastic 
scattering of 662 kev gamma rays", Proc. Phys. Soc, vol. 85, no. 6, 
pp. 1075-1081, 1965. 

[121] S. S. Nandi, S. K. Ghose, S. K. Sen Gupta, and N. Chaudhuri, "Study 
of anomalous dispersion in elastic scattering of 59.5 keV photons at 
K-absorption edges of target atoms", /. Phys. B, vol. 22, no. 8, pp. 
1175-1182, 1989. 

[122] A. Nath and A. M. Ghose, "Measurement of coherent scattering cross 
sections of gamma rays at small angles", Nucl. Phys., vol. 57, pp. 
547-564, 1964. 

[123] N. G. Nayak, K. M. Balakrishna, and K. Siddappa, "Experimental study 
on coherent scattering of 145.4 and 279.2 keV gamma rays by medium 
and heavy elements", J. Phys. B, vol. 26, no. 22, pp. 4117-4127, 1993. 

[124] N. G. Nayak, K. Siddappa, K. M. Balakrishna, and N. Lingappa, 
"Coherent scattering of 59.5 keV 7 rays by some medium and heavy 
elements", Phys. Rev. A, vol. 45, no. 7, pp. 4490-4493, 1992. 

[125] M. S. Prasad, G. K. Raju, K. N. Murty, V. A. N. Murty, and V. Lahshmi- 
narayana, "Elastic-scatering of 279 keV gamma-rays", Indian J. Pure 
Appl. Phys., vol. 16, no. 9, pp. 836-840, 1978. 

[126] M. S. Prasad, G. K. Raju, K. N. Murty, V. A. N. Murty, and 
V. Lakshminarayana, "Elastic scattering of 145 keV gamma rays", /. 
Phys. B, vol. 11, no. 23, pp. 3969-3976, 1978. 

[127] S. Puri, B. Chand, D. Mehta, M. L. Garg, N. Singh, and P. N. Trehan, 
"Differential cross section measurements for the elastic scattering of 
59.5 keV photons by elements in the atomic region 13 < Z < 82", 
Nucl. Instr. Meth. B, vol. Ill, no. 3-4, pp. 209-214, 1996. 

[128] R. Quivy, "Mesure absolue de la section efficace differentielle de 
diffusion coherente des gamma de 662 keV dans le plomb", Nucl. 
Phys. A, vol. 101, no. 2, pp. 417-423, 1967. 

[129] N. Ramanathan, T. J. Kennett, and W. V. Prestwich, "An experimental 
investigation of small angle photon elastic scattering", Can. J. Phys., 
vol. 57, no. 3, pp. 343-352, 1979. 

[130] D. V. Rao, "Measurements of whole-atom coherent and incoherent 
scattering cross-sections in some rare earth elements at low photon 
energies", Appl. Radiat. Isot., vol. 42, no. 9, pp. 855-859, 1991. 

[131] D. V. Rao, R. Cesareo, A. Brunetti, and G. E. Gigante, "Influence 
of solid state environment effects on measured elastic scattering cross 
sections in the X-ray regime and the associated anomalous dispersion", 
Phys. Scr., vol. 62, no. 1, pp. 81-87, 2000. 

[132] D. V. Rao, R. Cesareo, and G. E. Gigante, "Coherent and incoherent 
scattering of 14.93, 17.44 and 21.12 keV photons from Al, Cu, Sr, Cd, 
Ce, Pr, Sm, Pt, Au and Pb", Phys. Scr., vol. 50, no. 3, pp. 314-320, 
1994. 

[133] D. V. Rao, R. Cesareo, and G. E. Gigante, "Elastic and Compton 
scattering cross sections in the atomic region 12 ^ Z ^ 82", Appl. 
Radiat. Isot., vol. 47, no. 2, pp. 219-227, 1996. 

[134] D. V. Rao, R. Cesareo, and G. E. Gigante, "Rayleigh and Compton 
scattering cross sections for low, medium and high Z elements in the 
energy region 23.18 ^ E ^ 30.85 keV", Phys. Scr, vol. 54, no. 4, pp. 
362-367, 1996. 

[135] D. V. Rao, R. Cesareo, and G. E. Gigante, "Elastic and Compton 
scattering of X-ray photons from platinum in the momentum transfer 
region 1.3175 h A -1 < q < 2.0448fr A -1 ", Phys. Scr, vol. 55, no. 
3, pp. 305-309, 1997. 

[136] D. V. Rao, R. Cesareo, and G. E. Gigante, "Elastic scattering and the 
associated anomalous dispersion in the vicinity of Cd, In, Sn and Sb, 
K-absorption edges", Appl. Radiat. Isot., vol. 48, no. 6, pp. 789-806, 
1997. 

[137] D. V. Rao, R. Cesareo, and G. E. Gigante, "A new method to estimate 
elastic scattering cross sections in the X-ray region and the associated 
anomalous dispersion", Appl. Radiat. Isot., vol. 49, no. 7, pp. 835-844, 
1998. 

[138] D. V. Rao, R. Cesareo, and G. E. Gigante, "Elastic scattering and the 
associated anomalous dispersion in the energy range 8.63 ^ E ^ 42.75 
keV from heavy atoms", X-Ray Spectrom., vol. 27, no. 6, pp. 381-389, 
1998. 

[139] D. V. Rao and G. E. Gigante, "Coherent and incoherent scattering of 
low-energy X-ray photons in the atomic region 13^Z^82", Can. J. 
Phys., vol. 74, no. 1-2, pp. 10-16, 1996. 

[140] D. V. Rao, G. E. Gigante, and R. Cesareo, "Coherent and incoherent 
scattering of 42.75 and 47.24 keV x-ray photons scattered from Al, 



141 

142 
143; 
144 

145 

146 

147 

148 

149 
150 

151 
152 

153 

154; 

155 
156 

157; 

158; 

159 

160 

161 
162 
163 
164 

165 
166 
167 



Cu, Y, Mo, Au and Pb", X-Ray Spectrom., vol. 24, no. 4, pp. 172-176, 
1995. 

S.C. Roy, A. Nath, and A. M. Ghose, "Measurement of coherent scat- 
tering cross-sections of gamma rays at very low momentum transfer", 
Nucl. Instr. Meth., vol. 131, pp. 163-166, 1975. 

H. Schopper, "Die elastische streuung von 7-Strahlen bei kleinen 
streuwinkeln", Z. Phys. A, vol. 147, pp. 253-260, 1957. 

M. Schumacher, "Elastic scattering of 145-, 279-, 412-, and 662-keV 
7 rays from lead", Phys. Rev., vol. 182, no. 1, pp. 7-14, 1969. 
M. Schumacher, F. Smend, and I. Borchert, "Elastic scattering of 7- 
rays from lead for energies from 0.145 to 1.33 MeV", Nucl. Phys. A, 
vol. 206, no. 3, pp. 531-544, 1973. 

M. Schumacher and A. Stoffregen, "Rayleigh scattering of 59.54 keV 
photons from Zn, Mo, Sn, Ta, Au and Pb through angles from 60° to 
150°", Z. Phys. A, vol. 283, pp. 15-19, 1977. 

S. K. Sen Gupta, N. C. Paul, Bose J., G. C. Goswami, Das S. C, and 
N. Chaudhuri, "Atomic Rayleigh scattering of phonons in the vicinity 
of K-absorption edges", J. Phys. B, vol. 15, no. 4, pp. 595-602, 1982. 
S. K. Sen Gupta, N. C. Paul, S. C. Roy, and N. Chaudhuri, "New 
measurements of coherent and incoherent photon scattering by low, 
medium and high Z atoms", /. Phys. B, vol. 12, no. 7, pp. 1211-1224, 
1979. 

J. S. Shahi, S. Puri, D. Mehta, M. L. Garg, N. Singh, and P. N. Trehan, 
"Elastic scattering of 22.1 keV photons by elements in the atomic 
region 12 ^ Z ^ 92", Phys. Rev. A, vol. 55, no. 5, pp. 3557-3565, 
1997. 

J. S. Shahi, S. Puri, D. Mehta, M. L. Garg, N. Singh, and P. N. Trehan, 
"Large-angle elastic scattering of 59.54 keV photons by elements with 
12 ^ Z ^ 92", Phys. Rev. A, vol. 57, no. 6, pp. 4327-4334, 1998. 
K. Siddappa, N. G. Nayak, K. M. Balakrishna, and N. Lingappa, 
"Experimental studies on atomic form factors at 4.808-A -1 photon 
momentum transfer", Phys. Rev. A, vol. 39, no. 10, pp. 5106-5110, 
1989. 

F. Smend and H. Czerwinski, "Large-angle elastic scattering of 59.54 
keV photons by Kr and Xe", Z. Phys. D, vol. 1, pp. 139-140, 1986. 

F. Smend, M. Schumacher, and I. Borchert, "The Z-dependence of the 
elastic scattering of 7-rays", Nucl. Phys. A, vol. 213, no. 2, pp. 309- 
316, 1973. 

K. G. Standing and J. V. Jovanovich, "The elastic scattering of Co 60 
7-rays", Can. J. Phys., vol. 40, no. 5, pp. 622-653, 1962. 
R. B. Taylor, P. Teansomprasong, and I. B. Whittingham, "Elastic 
scattering of 152 Eu 7 rays by tungsten", Aust. J. Phys., vol. 34, pp. 
125-133, 1981. 

R. B. Taylor, P. Teansomprasong, and I. B. Whittingham, "Elastic 
scattering of 152 Eu and 154 Eu 7-rays", Nucl. Instr. Meth. A, vol. 255, 
no. 1-2, pp. 68-74, 1987. 

K. M. Varier and M. P. Unnikrishnan, "Elastic and inelastic scattering 
of 59.54 keV gamma rays", Nucl. Instr. Meth. A, vol. 280, no. 2-3, pp. 
428-432, 1989. 

W. M. J. Veigele, "Photon cross sections from 0.1 keV to 1 MeV for 
elements Z = 1 to Z= 94", Atom. Data Nucl. Data Tables, vol. 5, no. 

I, pp. 51-111, 1973. 

R. R. Wilson, "Scattering of 1.33 MeV gamma-rays by an electric 
field", Phys. Rev., vol. 90, no. 4, pp. 720-721, 1953. 
C. Gowda, T. K. Umesh, and R. Gowda, "A simple method to measure 
the total coherent scattering cross sections of elements", Nucl. Instr. 
Meth. B, vol. 95, pp. 127-130, 1995. 

C. Gowda and R. Gowda, "Studies on coherent scattering be- 
low 1 MeV", Thesis, Univ. of Mysore, 1994, Online. Available: 
http://dspace.vidyanidhi.org.in: 8080/dspace/handle/2009/l 328 
IEEE Computer Society, "IEEE Standard for Software Verification and 
Validation", IEEE Std 1012-2004, Jun. 2005. 

ISO/IEC, "International Standard, Information Technology Software 
Life Cycle Process, ISO 12207", 2008-2011. 

G. A. P. Cirrone et al., "A Goodness-of-Fit Statistical Toolkit", IEEE 
Trans. Nucl. Sci., vol. 51, no. 5, pp. 2056-2063, 2004. 

B. Mascialino, A. Pfeiffer, M. G. Pia, A. Ribon, and P. Viarengo, "New 
developments of the Goodness-of-Fit Statistical Toolkit", IEEE Trans. 
Nucl. Sci., vol. 53, no. 6, pp. 3834-3841, 2006. 
R. K. Bock and W. Krischer, "The Data Analysis BriefBook ", Ed. 
Springer, Berlin, 1998. 

A. N. Kolmogorov, "Sulla determinazione empirica di una legge di 
distribuzione", Gior. 1st. Ital. Attuari, vol. 4, pp. 83-91, 1933. 
N. V. Smirnov, "On the estimation of the discrepancy between empir- 
ical curves of distributions for two independent samples", Bull. Math. 
Univ. Moscou, 1939. 



31 



[168] T. W. Anderson and D. A. Darling, "Asymptotic theory of certain 

goodness of fit criteria based on stochastic processes", Anls. Ma. St., 

vol. 23, pp. 193-212, 1952. 
[169] T. W. Anderson and D. A. Darling, "A test of goodness of fit", JASA, 

vol. 49, pp. 765-769, 1954. 
[170] H. Cramer, "On the composition of elementary errors. Second paper: 

statistical applications", Skand. Aktuarietidskr. , vol. 11, pp. 13-74, pp. 

141-180,1928. 

[171] R. von Mises, "Wahrscheinliehkeitsrechnung und ihre Anwendung in 
der Statistik und theoretischen Physik", Leipzig: F. Duticke, 1931. 

[172] R. A. Fisher, "On the interpretation of x 2 from contingency tables, 
and the calculation of P", J. Royal Stat. Soc, vol. 85, no. 1, pp. 87-94, 
1922. 

[173] F. Yates, "Contingency table involving small numbers and the x 2 test", 

/. Royal Stat. Soc. Suppl, vol. 1, pp. 217-235, 1934. 
[174] K. Pearson, "On the x 2 test of Goodness of Fit", Biometrika, vol. 14, 

no. 1-2, pp. 186-191, 1922. 
[175] R. Barlow, "Statistics A guide to the use of statistical methods in the 

physical sciences", Ed: John Wiley & Sons, 1989. 
[176] S. C. Roy, L. Kissel and R. H. Pratt, "Elastic photon scattering at small 

momentum transfer and validity of form-factor theories", Phys. Rev. A, 

vol. 27, pp. 285-290, 1983. 
[177] T. V. George, L. Goldstein, L. Slama, and M. Yokoyama, "Scattering of 

Ruby-Laser Beam by Gases", Phys. Rev. 137, pp. A369-A380, 1965. 
[178] R. R. Rudder and D. R. Bach, "Rayleigh Scattering of Ruby-Laser 

Light by Neutral Gases" /. Opt. Soc. Am., vol. 58, pp. 1260-1266, 

1968. 

[179] S. J. Buckman, J. W. T. Cooper, M. T. Elford, M. Inokuti, Y. Itikawa, 

and H. Tawara, "Interactions of Photons and Electrons with Atoms", 

Springer- Verlag, 2000. 
[180] Shardanand and Y. Mikawa, "Photon Scattering Cross Sections at 

Lyman-Q! (1215.7 A) for He and Ne", /. Quant. Spectrosc. Radiat. 

Transfer, vol. 7, pp. 605-609, 1967. 
[181] Shardanand and A. D. P. Rao, "Absolute Rayleigh scattering cross 

sections of gases and freons of stratospheric interest in the visible and 

ultraviolet regions", NASA Technical Note TN D-8442, Washington 

DC, 1977. 

[182] P. Gill and D. W. O. Heddle, "Determination of the Refractive Indices 
of Gases in the Vacuum Ultraviolet. II The Rayleigh Scattering 
Method", /. Opt. Soc. Am., vol. 53, pp. 847-851, 1963. 

[183] P. D. Chopra and D. W. O. Heddle, "Polarization free measurements of 
Rayleigh scattering of Lyman ot\ J.Phys. B, vol. 7, no. 17, 2421-2428, 
1974. 

[184] J. P. Geindre, J. C. Gauthier and J. F. Delpech, "Rayleigh Light 
Scattering with a Low Power He-Ne Laser", Phys. Lett. A, vol. 44, 
pp. 149-150, 1973. 

[185] D. E. Gray, "American Institute of Physics Handbook", Third Edition, 

McGraw-Hill, pp. 6-110, 1972. 
[186] Allen, C. W, "Astrophysical Quantities" Third Edition, The Athlone 

Press, University of London, p. 92, 1973. 
[187] D. W. O. Heddle, R. E. Jennings and A. S. L. Parsons, "Determination 

of the refractive indices of gases in the vacuum ultraviolet. I. The 

Cerenkov radiation method", /. Opt. Soc. Am., vol. 53, pp. 840-846, 

1963. 

[188] G. I. Chashchina, V. C. Gladushchak and E. Y. Shreider, "Determina- 
tion of the Refractive Index of Argon and the Oscillator Strength of 
Its Resonance Lines", Opt. Spectrosc, vol. 24, pp. 542-543,1968. 

[189] R. B. Cairns, F. F. Marmo and J. A. R. Samson, "Photon Scattering 
by Argon in the Vacuum Ultraviolet", J. Opt. Soc. Am., vol. 60, pp. 
211-213, 1970. 

[190] M. Sneep and W. Ubachs, "Direct measurement of the Rayleigh 

scattering cross section in various gases", /. Quant. Spectrosc. Radiat. 

Transfer, vol. 92, pp. 293-310, 2005. 
[191] M. Han, C. H. Kim, L. Moneta, M.G. Pia, H. Seo, "Physics data 

management tools: computational evolutions and benchmarks" in Proc. 

Joint In. Conf. on Supercomputing Nucl. Appl. and Monte Carlo, p. 

10269, 2010. 

[192] M. Han, M. G. Pia, H. Seo, L. Moneta, and C. H. Kim, "Physics data 
management tools for Monte Carlo transport: Computational evolutions 
and benchmarks", in IEEE Nucl. Sci. Symp. Conf. Rec, pp. 95-101, 
2010. 

[193] M. Han et al., "New data libraries and physics data management tools", 
/. Phys. Conf. Ser, vol. 331, p. 042010, 2011. 



